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ABSTRACT 



Context. The analysis of unresolved stellar populations demands evolutionary synthesis models with realistic physical ingredients 
and extended wavelength coverage. 

Aims. To obtain a quantitative description of the first CO bandhead at 2.3 yum, to allow stellar population models to provide improved 
predictions in this wavelength range. 

Methods. We have observed a new stellar library with a better coverage of the stellar atmospheric parameter space than preceding 
works. We have performed a detailed analysis of the robustness of previous CO index definitions with spectral resolution, wavelength 
calibration, signal-to-noise ratio, and flux calibration. 

Results. We define a new line-strength index for the first CO bandhead at 2.3 yum, Deo, better suited for stellar population studies 
than previous index definitions. We compute empirical fitting functions for the CO feature as a function of the stellar parameters (T cff , 
logg and [Fe/H]), showing a detailed quantitative metallicity dependence. 

Key words, atlases - stars: fundamental parameters - globular clusters: general - galaxies: stellar content 
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One of the most important challenges in modern astrophysics 
is the proper understanding of the stellar content of unresolved 
systems, such as extragalactic globular clusters and galax- 
ies in different environm ents . Since the pioneering wo rk of 
ICrampin&Hovie1(ll961l) and iTinslevI (fl972l 1 19781 1 19801) . this 
has been accomplished through the comparison of the photo- 
metric and spectroscopic data with so-called evolutionary stel- 
lar population synthesis models, which make use of theoretical 
isochrones and libraries of spectral energy distributions (SEDs), 
either t heoretical, empirical or mixed (for more recent models 
see e g. lVazdekis et alJ2003HBruzual & Charlotl2003LlMarastonl 
2005). The most powerful approach to achieve this goal is to 
compare a number of observed line-strengths indices with their 
model predictions, providing in this way constraints to the rel- 
evant physical properties of the systems, namely age, metal- 
licity, initial mass function (IMF), and the relative abundance 
of different chemical species. Since, obviously, the reliability 
of model predictions improves as more realistic physical in- 
gredients are included, an important effort has been devoted to 
improve the quality of the SED libraries. Theoretical libraries 
usually exhibit systematic discrepancies among t hemselves and 
when comp ared with observational data (e.g., iLeieune et al.l 
1 19971 [19981) . Although the alternative empirical libraries consti- 
tute a coarse grained, and usually incomplete (especially for non 
solar metallicities and non solar abundance ratios) sampling of 
the space of stellar atmospheric parameters, the use of empirical 



20021) can help to reduce these effects (e.g.. 



1994; Cenarro et al 
IWorthevll 1 994 1 Vazdekis et"alll2003l) . 

Up to date, most of the observational effort has been 
focused to obtain complete libraries in the optical range. 
However, a full understanding of the physical properties of in- 
tegrated stellar systems cannot be achieved ignoring other spec- 
tral windows. In this sense, the CO features in the K band 
have been used by many researchers to i nvestigate the stel- 
lar content of galaxies, including ellipticals (|Fro gel et al. 19751 




980; 



Mobasher & James 1996 
Mannucci et al l 1200 it 



2000; James & Mobasher 



Silvaetal. 2008; Davidge et al 



spirals dJames & Seigarl ll999riBendo & JosephlhoO^ 
comp act galaxies dDavidge & Jensen! 120071: iMieske & Kroupa 
2008), starbursts and activ e galactic nuclei dDovon et al 
199A iRidgwav et all Il994t IShier etap I1996L [Puxley et al 



1997 



Goldader eFd]ll997t IVanzi & Rieke 1997; Mavya 1997; 



Ivanov et alj|2000l iHill etalj Il999t iRiffel et al.ll2007i among 



others. These strong absorptions are the ba ndheads formed in 
the fi rst overtone (Av = +2) bands of CO (Klein mann & Haffl 
1986). Despite the common use of these spectral features for 
stellar population studies, a proper characterization of the CO 
bands with stellar atmospheric parameters is still lacking. For 
that reason, we present in this work an improved study of the 
infrared region around 2.3 fim, where the first bandhead of the 
strong CO absorptions appear. In particular, we have observed 
a new library of stars which clearly surpasses preceding works 
(see § 12. U in the coverage of the stellar atmospheric parame- 



ters. After a thorough analysis of previous index definitions that 
have been used to measure the first CO bandhead, we present a 
new index, Deo, which is well suited for stellar population stud- 
ies. This new index depends very little on spectral resolution (or 
velocity dispersion), less sensitive to uncertainties in radial ve- 
locities, and can be measured with poorer S/N ratios. 

In Section |2 we present the new stellar library, highlighting 
the improvements over previous libraries, the sample selection 
as well as an overview of the observations and the data reduc- 
tion. A detailed discussion of the Deo index definition is given 
in Section [3] This section also includes a comparative study of 
the robustness of the new index to relevant effects. The measure- 
ments of the Deo index for the stellar library, and their associated 
error estimates appear in Section|4] Section[5]describes the stel- 
lar atmospheric parameters used to compute the fitting functions, 
which are derived in Section[6] Finally, Appendix A includes the 
tables with all the Deo measurements for all the stars used for the 
fitting functions, as well as their stellar atmospheric parameters. 



2. The new stellar library 

2.1. Previous work 



Several authors have compiled, for differ ent purposes, spec- 
troscopic stellar libraries in th e K band dJohnson & Mendezl 



1970[ 


Kleinmann & Hall 1986; Lancon & Rocca-Volmerange 


1992; 


Ali et al.l 11995b lHanson et al.l 1 19961 1 Wallace & Hinkle 


19961 119971 Ramirez et al.1 1 19971 iForster Schreiberl 


2000; 


Lancon & Woodll2000t llvanov et al]|2004t lHanson et al.l 


2005; 



Ranadaetal 2007b. TableU summarizes the previous stellar li- 



braries in the K band, including the number of stars, spectral 
range, spectral resolution and spectral types of the stars in each 
library. Due to the high S/ N ratio of their spect r a, it is interesting 
to highlight the library of iKleinmann & Halll (|19 86) (hereafter 
KH86), which cont ains 26 stars, but only with solar abundances. 
llvanov et al.l ([2004) present a library with 218 stars, which are 
not flux calibrated. The poor metallicity coverage for these pre- 
vious libraries (see Figs. Q] and |2]l has not made it possible ex- 
plore the metallicity dependence of the spectral features in the K 
band. 

2.2. Sample selection 

We have observed a new stellar library in the K band which 
comprises 220 stars. The observed sample is a subset of 
MILES (Medium-re solution Isaac Newton Telescope Library of 
Empi rical Spectra; Sanchez-Blazq uez et al.ll2006l ICenarro et al.l 
2007), a stellar library in the optical r ange with well know n at- 
mospheric parameters for all the stars (Cen arro et al.l l2007). Our 
final stellar sample includes stars in the following stellar param- 
eter ranges: 

2485 K < r eff < 13404 K, 
-0.34 dex < logg < 5.30 dex, 
-2.63 dex < [Fe/H] < +0.98 dex, 

where [Fe/H] = logZ - logZ . 

This library clearly has a larger metallicity coverage than the 
previous ones (see TableQ]and Fig. |2]for a comparison between 
diffe rent works), and contains 8 stars in co mmon with KH86, 2 3 
with lWallace & Hinkld d 19971) and 39 with llvanov etail (l2004 . 



Table 2. Observational configurations. 



Telescope 


CAHA 3.50 m 


TNG 3.56 m 


Instrument 


Q-CASS 


NICS 


Slit width (") 


0.60 


0.75 


Grism 


#1 


KB 


Filter 


K 


— 


Spectral coverage 


2.01-2.43 yum 


1.95-2.34 yum 


Dispersion 


2.527 A/pix 


4.375 A/pix 


FHWM 


6.8 A 


11.3 A 


Detector 


Hawaii-I 


Hawaii-I 



2.3. Observations and data reduction 

The bulk of the stellar library (217 stars) was observed dur- 
ing a total of 13 nights in five observing runs from 2002 to 
2005 on the 3.5 m telescope at Calar Alto Observatory (CAHA, 
Almeria, Spain) with Q-CASS. A subsample of the stellar li- 
brary (52 stars) was observed again at the Telescopio Nazionale 
Galileo (TNG) at Roque de los Muchachos Observatory (La 
Palma, Spain) with NICS (Near Infrared Camera Spectrometer) 
in February 2006 and May 2007, plus 3 new stars. The details of 
the instrumental configuration for both runs are given in Table|2] 

Each star was observed several times at different positions 
of the slit (standard procedure for infrared observations) to per- 
form a reliable sky subtraction. Halogen lamps (on and off) and 
arc lamps were observed for flat-fielding, and C-distortion cor- 
rection and wavelength calibration, respectively. Vega type (AO) 
stars were observed at different airmasses during each night in 
order to calibrate in relative flux and eliminate telluric lines in 
the stellar spectra. 

We carrie d out a standa rd data reduction in the infrared us- 
ing REDmE dCardielll 19991) . a reduction package which allows 
a parallel treatment of data and error frames. The reduction pro- 
cess includes flat-fielding, sky subtraction by subtracting consec- 
utive images (A-B), cosmetic cleaning, C-distortion correction 
and wavelength calibration with arc lamps, spectrum extraction 
and relative flux calibration. Atmospheric extinction was cor- 
rected by using extinction coefficients (namely the relative con- 
tributions of the Raleigh scattering and the aerosol ext i nction ) 
derived for CAHA Observatory by iHopp & Fernandez! (|2002). 
Those coefficients were extrapolated for La Palma Observatory 
to correct the stars observed in this observatory. 

Some of the reduction steps that requires more careful work 
are explained in detail in the following subsections. 

2.3.1. Wavelength calibration 

Arc spectra of Argon lamps were acquired to perform the C- 
distortion correction and the wavelength calibration. Due to in- 
stabilities and flexures of the instrument and the telescope, cal- 
ibration arc frames were obtained after each star observed at 
CAHA. In the K band, the typical number of known lines in 
the arc spectrum is rather low (just six in our instrumental con- 
figuration). Because of that, the wavelength calibration was not 
accurately enough and a second order wavelength correction was 
performed b y identifying OH air-glow lines in the sky spec trum 
of each star dOliva & Origlialll992h iRousselot etakl feoOO). For 
observations at CAHA, a polynomial fit of the differences be- 
tween the observed and theoretical OH lines for the sky spec- 
trum was necessary. The wavelength calibration polynomial is 
expressed as a function of position as 

W(x) = A(x) + z(A(x)), (1) 



02 



Of) 
O 




[Fe/H] 

□ >0.25 
(-0.25. + 
(-0.75,-0 
(-1.25,-0 
it (-1.75,-1 
£> <-1.75 
* unknown 



25) 
25) 
75) 
25) 



0.5 



1 



1.5 



( = 5040/T eff ) 

Fig. 1. \ogg - 9 diagram for the stellar library presented in this work, where 6 = 5040/ T^. Different symbols are used to indicate 
stars of different metallicities, as shown in the key. The boxes display the regions of the corresponding local fitting functions for the 
new CO index (see § l6.2.2l for details). 
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Fig. 2. Stellar parameter coverage in the stellar libraries of KH86, IWallace & Hinkiel (Il997h . Ilvanov et al.1 (|2004), and this new 
library. We present separatel y dwarf and giant (and su pergiant) stars (up per and lower panels, respecti vely). The stellar parameters 
of the stars from KH86 and Walla ce & Hinkiel (\991) were taken for Ca vrel de Strobel et alj d200ll) if available. Otherwise, we 
assigned solar metallicity and T e ff from their spectral type following the tables of|Langl(ll99ir(sinall asterisks). 



Table 1. Main characteristics of previous spectroscopic stellar libraries in the K band and the new stellar library presented in this 
work. 



Reference 


Number of 


Spectral range 


Spectral resolution 


Spectral 


Notes 




stars 


w 


(R = A/AA) 


types 




Johnson & Mendez (1970) 


32 


1.2 


-2.5 


550 


A-M, I-V 


Low resolution 


Kleinmann & Hall (19861 


26 


2.0 


-2.5 


2500-3100 


F-M, I-V 


Solar abundances 


Laneon & Rocca-Volmeranee (1992) 


56 


1.4 


-2.5 


550 


O-M, I-V 


Low resolution 


Alietal. (1995) 


33 


2.0 


-2.4 


1380 


F-M, V 


Dwarf stars 


Hanson et al. (1996) 


180 


2.0 


-2.2 


800 - 3000 


O-B, I-V 


Hot stars, not CO region 


Wallace & Hinkle (1996) 


12 


2.02 


-2.41 


> 45000 


G-M, I-V 


Few stars 


Ramirez et al. (1997) 


43 


2.19 


-2.34 


1380, 4830 


K-M, III 


Giant stars 


Wallace & Hinkle (1997) 


115 


2.0 


-2.4 


3000 


O-M, I-V 


Solar abundances 


Forster Schreiber (2000) 


31 


1.90 


-2.45 


830, 2000 


G-M I-III 


Giant and supergiant stars 


Lancon & Wood (2000) 


77 


0.5 


-2.5 


1100 


K-M, I-III 


Giant and supergiant stars 


Ivanov et al. (2004) 


218 


1.48 


-2.45 


2000 - 3000 


G-M, I-V 


Not flux calibrated 


Hanson et al. (2005) 


37 


2.0 


-2.2 


8000-12000 


O-B, I-V 


Hot stars, not CO region 


Cushing et al. (2005) 


26 


0.6 


-4.1 


2000 


M-T, V 


Extremely cold dwarf stars 


Ranada et al. (2007) 


114 


2.05 


-2.19 


2200 


O-M, I-V 


Not CO region 


This work 


220 


2.11 


-2.37 


2500 


O-M, I-V 


Improved metallicity coverage 



being A (x) the initial approximation to the wavelength calibra- 
tion 



A (x) = 2^ a i x \ 



(2) 



where x is the position in the spectral direction, and a-, are the 
coefficients of the n th -degree calibration polynomial. The second 
order correction is given by 



z(A(x))=J]bjA j , 



(3) 



where z (A) are the computed differences between the observed 
and theoretical OH lines, and bj are the coefficients of a 
m th -degree polynomial. In our case, it is a second order poly- 
nomial. The final wavelength calibration polynomial correction 
is 



z(x) 



2> 

./=() 



(2> 



k=0 



x J \ 



(4) 



where q are the coefficients of the new correction polynomial of 
order I = m ■ n as a function of the position x. 

In the case of the observations at the TNG, we compared 
the observed sky spectrum with the well calibrated sky spectrum 
from CAHA and we observed constant wavelength differences 
between them. We cross-correlated both spectra and we applied 
this constant shift to the wavelength calibration of TNG spectra. 

Finally, we c hecked the final spectra w ith K band spectra 
from KH86 and IWallace & Hinklel (Il997l) . since they used a 
Fourier transform spectrometer, which implies a very accurate 
wavelength calibration of the spectra. Although in the case of ob- 
servations at CAHA no differences were obtained, for the TNG 
observations we had to apply a constant shift in order to achieve 
the correct wavelength calibration. The origin of this discrep- 
ancy is found in the lack of OH sky lines in the reddest wave- 
length region of the K band, which prevented from an accurate 
cross-correlation of the sky spectra in that region. 
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-0.5 




Fig. 3. Upper panel: Histogram with the different K factors (see 
Eq-Eh employed to correct a particular stellar spectrum from tel- 
luric absorptions. Lower panel: Example of a flux standard spec- 
trum before (dotted line) and after (solid line) telluric correction 
(see details in the text). 



2.3.2. Flux calibration and telluric correction 

There are two ways of flux calibrating infrared spectra. The first 
method consists of observing a solar type star close to the star 
to be calibrated. The solar type star is red uced as usual, dividing 
the fi nal spectrum by the solar spectrum dLivingston & Wallace! 
Il99ll) degraded to the same spectral resolution as the problem 
star. In this way, a spectrum with the information about the re- 
sponse curve and the telluric lines is obtained. This spectrum is 
then rectified by the ratio between the blackbody spectra at the 



temperature corresponding to the solar type star and the solar 
temperature in order to obtain the correct continuum. This final 
spectrum is used to (relative) flux calibrate and carry out the tel- 
luric correction in the stars to be calibrated. The advantage of 
this method is that it is easy to find a star of this type for each 
observation. 

A second method consists in the observation of Vega type 
stars at different airmasses during the night. The main reason 
for choosing these stars is that they are known to have no rele- 
vant features in our observational window, except the Bry line. 
After wavelength calibration, Vega type stars are divided by the 
well known theoretical Vega spectrum in our spectral range. In 
that way, a spectrum with both the response curve and the tel- 
luric corrections is obtained. The final stellar spectrum is then 
obtained after dividing each star by this spectrum. In this work, 
we have used this second approach. 

2.3.3. Second order telluric correction 

As we mentioned in the previous section, we used the flux stan- 
dard star in order, not only to flux calibrate the stellar spectra, but 
also to correct simultaneously for the telluric absorption lines. 
Due to the variability of the observing conditions during the 
night, some telluric lines are badly corrected by applying the re- 
sponse curve derived from the flux standard star. For that reason 
an extra correction was necessary. First of all, we computed a 
reference spectrum with the information of the telluric lines. For 
observations at CAHA, we checked the response curves for each 
night looking for the spectrum with the best removal of telluric 
features. The ratio between each response curve and the previ- 
ous spectrum free from telluric contamination provides the tel- 
luric spectra that we used to correct all the stellar spectra. In the 
case of the TNG observations, the telluric spectrum was obtained 
by dividing the flux standard spectra at high and low airmasses, 
observed during each night. The telluric spectrum in both obser- 
vatories, obtained as explained above, consists mainly of differ- 
ences in the strength of the telluric absorption lines. To correct 
for these lines, we modified their intensity by multiplying by an 
adjustable factor K, i.e., 

s* = (s -i)xjr + i, (5) 

where Sq is the telluric spectrum, and S k is the telluric spectrum 
adjusted to correct a specific stellar spectrum. We divided the 
latter stellar spectrum by S k and computed the r.m.s. in the cor- 
rected spectrum. The best correction factor, K, is the one which 
minimizes the r.m.s. This method was applied to different iden- 
tified telluric lines separately, since they do not vary in the same 
way. This effect can be seen in the histogram of Fig[3] top panel, 
where the number of telluric lines corrected by a factor K for 
a given spectrum is represented. In Fig. [3] bottom panel, we 
present an example of the telluric lines correction for a given 
flux standard spectrum. Notice that the telluric absorption lines 
can be present even after flux calibration and it is important to 
correct them in infrared spectroscopy. 

3. Index definitions for the CO band at 2.3 fim 

3. 1 . The K band region 

The most prominent features in the K band are due to the 
rotational-vibrational transitions of the CO molecule around 
2.3 fim. Important absorptions are also produced by other metal- 
lic species, such as Na I, Fe I, Ca I and Mg I (KH86, see TableO. 



Tabl e 3. Main spectros copic features in the K band (from 
lKleinmann&Halllll986l) . 



Species A Transition Lower state 



(fim) energy (eV) 



HI Bry 


2.1661 


4 I F l > - 7 2 G 


12.70 


Nal 


2.2062 


4 S 2 S l/2 -4p 2 P° l/2 


3.19 


Nal 


2.2090 


4 S 2 S l/2 -4p 2 P° /2 


3.19 


Fel 


2.2263 


x 5 F° - e 5 D 3 


5.07 


Fel 


2.2387 


x 5 F* - e 5 D 2 


5.04 


Cal 


2.2614 


4d 3 D 3XX - 4fF° 4 


4.68 


Cal 


2.2631 


4d 3 D 32i -4fF° 


4.68 


Cal 


2.2657 


4d 3 D 3Xl -4fFl 


4.68 


Mgl 


2.2814 


4d 3 D 3Xl - 6fF° 34 


6.72 


12 CO(2,0) 


2.2935 


(2,0) bandhead' ' 


0.62 


12 CO(3,l) 


2.3226 


(3,1) bandhead 


0.86 


13 CO(2,0) 


2.3448 


(2,0) bandhead 


0.32 


12 CO(4,2) 


2.3524 


(4,2) bandhead 


1.12 



The only hydrogen line in this spectral range, Bry, is present 
generally in absorption in O and B stars, going into emission for 
hig h-luminosity stars ( for a detailed study of this kind of stars, 
see lHanson etatll 19961) . 

In this paper, we focus our study on the first CO bandhead 
at 2.29 fim. Contrary to the Cal and Nal, the contribution of 
other species to the CO absorption is almo st negligible (see 
I Wallace & Hinkld 1 19961: iRamirez et all 1 19971 for a further dis- 
cussion). 

3.2. Previous definitions 

In order to measure the CO absorption at 2.3 fim in an objective 
way, several authors have proposed different index definitions. 
Bald win et alj d 19731) suggested a photometric system to mea- 
sure the CO features based on two narrow filters (AA = 0. 10/im) 
centered at 2.30 fim and at 2.20 //m for the CO absorption and the 
continuum, respectively. The CO index was defined as the differ- 
ence of the two filters relative to the values obtained for a Lyrae, 
in magnitudes. Following this idea, Frogel et al d 19781) defined 
the most used photometric CO index (CO p hotX with slightly dif- 
ferent filter parameters (AA = 0.08 /urn, for the CO filter centered 
at 2.36 fim, and AA — 0.1 1 fim for the filter centered at 2.20 fim 
for the continuum estimate). 

The first spectroscopic CO index (CO^ 8 ) for the CO(2,0) 
bandhead at 2.3 fim was defined by KH86 as 

CO™ g = -2.5 log COkh = -2.5 log p, (6) 

where COkh = Ta/Fc is the ratio between the fluxes in- 
tegrated over narrow wavelength ranges centered in the ab- 
sorption line (/LI2.29305 - 2.29832 fim) and the nearby con- 
tinuum (AA2. 28728 - 2.29252 fim), measured in magnitudes. 
These band limits ha ve been used to me asure the index as an 
equivalent width (e.g. Origl ia et al.lll993[). B oth measurements 
can be converted using dOriglia & Olivall2000l) 

rn mag /, W, (2.29) \ 

where CO^ 8 is the spectroscopic index initially defined by 
KH86 measured in magnitudes, and W,j (2.29) is the same in- 
dex measured as an equivalent width. 



Table 4. Spectroscopic CO index definitions. References are 



CO™ 8 



4. Sp> 
; dKleii 
(Froge 



Kleinmann & Halll 19861) . Ipu X iey dPuxley et al.ll997h and 



iFiogei dFrogel e t al. 2001). Wavelengths are in vacuum. 



Index 



-co^ 

Ipuxley 



Continuum 
bands (/jm) 



Absorption 
bands (/Jin) 



Comments 



Frogel 



2.2873-2.2925 
2.2530-2.2610 
2.2700-2.2780 
2.2850-2.2910 
2.2300-2.2370 
2.2420-2.2580 
2.2680-2.2790 
2.2840-2.2910 
2.2460-2.2550 
2.2710-2.2770 



2.2931-2.2983 
2.2931-2.3200 



Color-like index 
Generic index 



2.2910-2.3020 Generic index 



2.2880-2.3010 



Generic 
discontinuity 



iDoyon et al.l ([1994) studied the behaviour of CO p hot and in- 
dicated several reasons to introduce their new spectroscopic def- 
inition 



CO sp = -2.5 log < fl 2 .36 >, 



(8) 



where < R236 > is the mean value of the rectified spectrum (nor- 
malized in the continuum) in the 2.31 - 2.40 fim range. This 
rectified spectrum is obtained by fitting the continuum in the 
2.00 - 2.29 /mi range with a power law (F^ oc A~ a ), due to the 
similarity of t he stellar spectrum in the K band to a Rayleigh- 
Jeans law. As lOriglia & Olival d2000l) indicated, this index is just 
the equivalent width over the CO range relative to a continuum 
which is extrapolated from shorter wavelengths. As a main ad- 
vantage, this definition allows measurements of the CO even 
from poor quality spectra. 

Ot her authors (for example. lLancon & Rocca-Vormera nge 

Il992t lRamirezetailll997t iForster Schreiberl l2000h proposed 
their own definitions, adopting the bandpasses for the absorp- 
tion and the continuum without conside ring the use of thos e def- 
initions in general situations. Recently, Rif fel et ail J2007) mea- 
sured the CO absorption at 2.3 /mi as an equivalent width be- 
tween /i/12.2860 - 2.3100 /mi, computing a continuum defined 
as a spline using points free from emission/absorption lines in 
the broad interval i^2.2350 - 2.3690 /mi. 

After a study of different band limits for the CO index mea- 
sured in terms of equ i valent widths for stars of different spec- 
tral type, Puxl ey et al.1 {l997) proposed to extend the absorption 
band of KH86 up to the end of the CO (2,0) band (2.320 /im), 
and the use of three different bands to estimate the continuum 
(AA2.253 - 2.261 /mi, 2.270 - 2.278 //m and 2.28 5 - 2.291 urn) 
Note that this type of index definition is what Cenarro et al.l 



(2001a) called a generic index. iPuxlev et all d 19971) adopted this 
definition because it allows to distinguish between giant and 
supergiant stars and the correction for velocity dispersion is 
smaller than with the other definitions. 

An additional definition was introduced by IFrogel et alj 
(2001), in which the CO absorption feature is measured using 
multiple bandpasses to estimate the pseudo-continuum level. 

3.3. New index definition 

Even though the number of different CO index definitions is 
large, we have explored in detail whether any of these is actually 
well suited for a practical study of this spectroscopic feature in 
the integrated spectra of galaxies. Curiously, from t he list of pre- 
vious index definitions only the one presented by IPuxlev et al.l 




wavelength (/im) 

Fig. 4. Limits of the bandpasses in the definitions for the 
CO index propose d by: a) KH86; b) IPuxlev et all dl997l) : c) 
IFrogel et all d200ll) : and d) the new CO index presented in this 
work. Grey and open bands represent absorption and contin- 
uum bandpasses, respectively, for each index definition, super- 
imposed on the spectrum of HD 137704. 



d 19971) (based on the previous definition by KH86) was designed 
taking into account the variations of the index with radial veloc- 
ity and velocity dispersion, both of them very important in the 
study of galaxies. In this work we have spent an additional ef- 
fort to investigate the possibility of finding an optimal CO index 
definition that could improve all the previous definitions. 

In order to carry out this task, we have focused our efforts on 
defining a CO index which is less sensitive to low signal-to-noise 
ratios, degradation due to spectral resolution and/or velocity dis- 
persion, and errors in wavelength calibration (or errors in radial 
velocity), and in relative flux calibration (see the next sections 
for a further study of each case). 

After exploring different possibilities for the definition of 
the new index, we propose to measure the CO at 2.29 /mi as a 
generic discontinuity, i.e., as the ratio between the average fluxes 
in the continuum and in the absorption bands 

£ ' F cJ (A) dA 

, = 1 JAci, 



/ , (^e,/ 2 - A c j { ) 



D 



generic — 



§1 



(9) 



F a AX) dA 



/ t (A1J2 - ^-a,h) 



where D gener i c is the generic discontinuity, and F a j (A) and 
F c j (A) are the flux in the n a absorption bands and n c continuum 
bands, respectively. Finally, A x j t and A x j 2 are the lower and upper 
wavelength limits of the / th band x (where x is a or c). This new 
definit ion is similar to the B4000 index defined by Gorga seit al.l 
(1 19991) but using more than one bandpass to define the contin- 
uum and the absorption regions. 

Here we propose to measure the CO feature at 2.3/mi as a 
generic discontinuity, Deo, using two bandpasses for the con- 
tinuum (n c = 2) and one bandpass for the absorption region 



(n a = 1). The limits of these bands are also listed in TableHJ We 
selected the number of bandpasses and their location taking into 
account several factors. Concerning the continuum bandpasses, 
we have eluded the Ca I and Mg I features, trying not to extend 
too far towards shorter wavelengths in order to avoid potential 
systematic effects arising in the flux calibration of wide line- 
strength indices. In the case of the absorption region, one single 
bandpass is enough to cover the first CO bandhead. Compared 
with previous definitions, we decided to shift slightly the blue 
continuum bandpass limit to obtain an index which is more sta- 
ble with velocity dispersion (i.e. spectral resolution) and with 

radial velocit y uncertainties . 

Following Cardi el et al1 (fl998). it is not difficult to show that 
the expected variance in a D genel i c index can be computed as 

<r 2 [D generic ] = ^— (10) 

'a 

where is the total flux per wavelength unit in the continuum 
(x = c) and the absorption (x = a) region, determined from the 
coaddition of the flux in all the corresponding bandpasses, i.e., 

H x y pixels 

r x = &^-^ , (ii) 

i=\ 

being n x the number of bandpasses in either the continuum or the 
absorption region, is the linear dispersion (in A/pixel), N 1 .^ 

is the number of pixels covered by the / th bandpass of the x re- 
gion (with x equal to c or a), and Ak is the central wavelength of 
the k th pixel. The variance in these total fluxes are simply com- 
puted as the quadratic sum of the individual variances in each 
pixel, i.e., 

N' 

n x pixels 

crl = &±LJz± (12) 

* r n x ] z 

. (=i 

where, in particular, cr 2 F (A*) is the variance corresponding to 

the random error in the k th pixel. It is important to highlight that 
in these expressions we are assuming that the random errors in 
each pixel are not correlated. 

3.4. Sensitivities of the indices to different effects 

In this section, we discuss the sensitivity of previous spectro- 
scopic indic es defined by KH86 (CO™ 8 ), iPuxlev et all d 19971) 
(Ipuxiey) and iFrogel et alJ (120011) (IFrogel), and the new CO in- 
dex (Deo) to velocity dispersion (or spectral resolution), wave- 
length calibration (radial velocity), relative flux calibration and 
signal-to-noise (S/N) ratio. A fifth index, Dprogei, is considered: 
a generic discontin uity based on the same bands proposed by 
Fro g el et"ai] (|2001). In Fig. [4] we show the bandpasses for these 
index definitions. For this study, we select ed from the high res- 
olution library of Wallace & Hinkle (1996) three stars with sim- 
ilar spectral type (M2-5; chosen because of their strong CO 
features) and different luminosity class (supergiant, giant and 



dwarf) in order to account for differences depending on the type 
of star. The resolution of the spectra is 0.54 A(FHWM) and they 
are all shifted to rest-frame. 

3.4.1. Spectral resolution and velocity dispersion broadening 

In order to study the sensitivity of the spectroscopic CO indices 
to the spectral resolution or velocity dispersion broadening (cr), 
we broadened the stellar spectra of the selected stars with addi- 
tional cr's from the initial <x () up to cr = 400 km/s (in steps of 10 
km/s). The different indices were measured on all these broad- 
ened spectra and we computed the ratio between the index (I) 
at each cr and the index measured on the original spectrum (Io). 
Fig.[5](left column) shows this ratio as a function of the velocity 
dispersion for the definitions we are studying. Compared to the 
other index definitions, the two generic discontinuities (Dprogei 
and Deo) are clearly the less sensitive to velocity dispersion 
broadening. 

3.4.2. Wavelength calibration 

Sometimes, errors in the wavelength calibration arise in the 
spectra even after a very careful reduction or due to an inaccu- 
rate radial velocity (v r ) estimate of the studied object. Because of 
that, it is important to define indices with the least possible sen- 
sitivity to this kind of uncertainties. To quantify this effect, we 
measured the CO absorption with the different index definitions 
in the stellar spectra shifted from -200 to +200 km/s with steps 
of 4 km/s in radial velocity. In Fig.[5](central column) we present 
the ratio between the index, I, measured at v r and the initial value 
Io (assumed v, = km/s) as a function of the considered v r for 
different types of stars. It is apparent from the figures that the 
indices CO^ 8 , Ip ux iey and IFrogel are very sensitive to radial ve- 
locity uncertainties, while Dp roge i and the new index definition 
Deo are more robust to this effect. 

3.4.3. Flux calibration 

As we explained in § 12.3.21 it is common to use theoretical spec- 
tra to recover the real shape of the continuum. This practise im- 
plies the knowledge of the temperatures of the standard stars. 
For that reason, we have studied the impact, during flux cali- 
bration, of an error in the temperature estimate of the standard 
stars. To analyse the impact when solar-type stars are used as flux 
standards, we computed the blackbody spectrum in the interval 
5600 < r e ff ^ 6300 K, and derived the ratio between these spec- 
tra and the blackbody, at solar temperature. To study the effect 
when Vega type stars are used as calibrators, we analysed the dif- 
ferences from the theoretical spectrum of Vega (T e s ~ 9400 K) 
and the real temperature of the Vega type stars (from 8400 to 
14400 K for our study). In both cases, we found that the changes 
in the continuum produced by differences in the assumed tem- 
perature of standard stars produce negligible differences in the 
measured indices. 

Finally, we studied the impact of a wrong curvature in the 
response curve, which is a typical source of systematic error. In 
order to obtain an estimate of this effect, as a first order approach 
we have artificially modified the continuum shape of the original 
spectra by multiplying them by a second order polynomial. This 
polynomial was chosen to pass through 3 fixed points, two at 
the borders of the wavelength range (where the polynomial were 
forced to be equal to 1 .0), and another point at the center of that 
range (where the polynomial was set to a variable parameter f3 
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Fig. 5. Study of the sensitivity of the different CO index definitions (from top to bottom, CO^ s , Ip ux iey> iFrogeb DFrogel and Dco) to 
several relevant parameters. Left column: ratio between the index I measured on the spectra broadened wit h a given cr, and the inde x 
In measured on the initial spectra (crn = 3 km/s, corresponding to the resolution of the stellar library of Wallace & Hinklel [T996). 
Central column: ratio between the index I for a given v r and the index In in the original spectrum (v r = km/s). Right column: ratio 
between the index I measured on the spectrum multiplied by a curved spectrum parametrized by /3 (see § l3.4.3l for details) and the 
index measured over the original spectrum, In, as a function of the curvature parameter /3. 



ranging from 0.5 to 1.5). In Fig. [6] we show different examples 
of these polynomials for distinct values of p. Note that, with this 
exercise, we simply studied the effect of a low frequency error 
in flux calibration. 

In Fig|5](right column) we present the ratio between the mea- 
sured index in the stellar spectrum multiplied by the polynomial 
for a given value of I, and the original one, i.e., In for/? =1.0 
(no additional curvature), as a function of the parameter /3. The 
sensitivity of each index definition to the (3 parameter is, not sur- 
prisingly, dependent on the location and wavelength coverage of 
the index bandpasses, and also depends on the way the pseudo- 
continuum is determined, and the absolute value of the index. 
For these reasons, CO^ E , Ip ux iey and lF r0 gei are the most sensi- 
tive to an error in the response curve. In particular, the CO^ 8 
depends strongly on the strength of the CO absorption. On the 
other hand, Ip ux i e y and lF ro gei> defined as generic indices, extrap- 
olate the continuum value into the absorption band and, in that 
way, the wrong curvature too. In addition, the generic discon- 
tinuities Dpiogei and Dco, computed as the averaged flux in the 



continuum and absorption bands, exhibit no differences between 
dwarf, giant and supergiant stars. 

3.4.4. S/N ratio 

One important issue to take into account in the definition of a 
new index is the dependence of the relative error of the measure- 
ments on the S/N ratio. In this sense, the aim is to find an index 
definition which provides the lowest relative error in the mea- 
surements with the lowest signal-to-noise ratio in the spectra. 
For that reason, we have studied the behaviour of relative errors 
measured with previously analyzed CO index definitions as a 
function of the S/N ratio. Using a particular stellar spectrum, we 
have simulated a set of one hundred spectra (and their associated 
error spectra) with random S/N(A) ratios in the ra nge 1.0-100.0 . 
For this task we have used the program index iFl (lCardiell2007l) . 
In Fig [7] we compare the results obtained for a giant star (same 



1 http://www.ucm.es/info/Astrof/software/indexf/ 




2.1 2.15 2.2 2.25 2.3 

A (yum) 

Fig. 6. Examples of the curved spectra used in the study of the 
sensitivity of the different index definition to a wrong curvature 
of the spectrum. The parameter /? determines the distance to the 
maximum/minimum in the center of the spectrum (shown here 
for the polynomial displayed with a solid line). 
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Fig. 7. Example of the study of the dependence of the relative er- 
ror in the measure of the first CO bandhead using different index 
definitions at a given S/N ratio per A. We represent the relative 
errors vs. S/N ratio (in logarithmic scale) measured on the sim- 
ulated spectra of a giant star (the results are independent of the 
type of star). Different symbols are used to indicate different in- 
dex definitions, as shown in the key. See text for details. 



results are obtained for supergiant and dwarf stars). Not surpris- 
ingly, the relative errors in all the definitions follow 



«r = 



S/N(A) 



(13) 



where c is a constant that depends on the particular index. 
This result was already found for atomic and molecular in- 
dices ( Cardie 1 et aTT l998). and for generic indices dCenarro et alj 
2001a). It is clear from Fig. [7]that the same holds for generic dis- 
continuities. Considering Eq.[l3] it is evident that, at a given S/N 
ratio, the lower relative errors correspond to the index definitions 
with lower c values. Table [5] list these values for the five index 
definitions under study. From these numbers and the data dis- 
played in Fig. [7] it is obvious that Deo is comparable to Dprogeb 
while CO'j^ 8 , Ipuxiey and iFrogei provide larger relative errors for 
a given S/N ratio. 



Table 5. Values of the constant c of Eq.[l3]for the different index 
definitions analyzed in § 13.41 



Index 


c 


C0 KH 


0.7537 


Ipuxiey 


2.0258 


Iprogel 


0.8123 


Dprogel 


0.1075 


D ro 


0.1198 



3.4.5. The best index definition 

Once we have studied the behaviour of the different CO index 
definitions as a function of all the relevant parameters, we can 
conclude that the Deo index definition is, in general, preferable. 

On one hand, CO^ s , I^gei and Ip ux iey are too sensitive of 
spectral resolution and errors in wavelength calibration and ra- 



and 



dial velocity. In addition, the behaviour of CO KH ", IFrogei 
Ipuxiey are also too sensitive to uncertainties in the spectropho- 
tometric system (i.e., flux calibration). 

When the sensitivity to S/N is included in the comparison, it 
is clear that the best definitions are the two generic discontinu- 
ities, namely Dprogei and Deo- Since the use of generic discon- 
tinuities for measuring the CO absorption is introduced for the 
first time in this paper, and considering that the Deo is practically 
insensitive to spectral resolution (or velocity dispersion broaden- 
ing) up to <x ~ 400 km/s, we propose to use the new definition, 
especially for the future analysis of integrated spectra. 

3.5. Conversion between different CO index systems 

In this section we give the calibrations to convert between the 
new CO index defin itio n and the CO indice s defined by KH86, 
iPuxley et all (1 19971) and lFrogel et alJ (120011) . In order to obtain 
these conversions, we have measured the indices on the sub- 
sample of stars observed at the TNG (3200 < r eff < 9625 K, 
0.00 < logg < 5.00, -1.73 < [Fe/H] < +0.36). The cali- 
brations were computed by deriving a least squares fit to the 
data. The fits are comp letely compatible with i ndex measure- 
ments on the KH86 and Wall ace & Hinkld (1 19971) spectra which 
are on the same spectrop hotometric system. Just six stars from 
IWallace & Hinkld d 19971) deviate more than 3cr from the fitted 
relation due to a problems in the continuum and the telluric cor- 
rection of those spectra. In Fig. [8] we show all these fits and the 
data used to compute them. 

The conversion between the index defined by KH86, CO^ g , 
and the new index Deo is given by 



Deo = 1.0407 (±0.0021) + 0.2317 (+0.0035) CO™ g 



(14) 



with r 2 = 0.9863. 

The expression to compute Deo from the index defined by 
Fro gel et alj d2001l) is 



Deo = 1.0507 (±0.0031) + 0.0077 (±0.0005) I Fro gei + 
+0.00007 (±0.00002) I^ rogel 



(15) 



with r 2 = 0.9802, where lF loge i is measured as an equivalent 
width (A). 

As we mentioned before, Ip ux iey is also measured as an equiv- 
alent width. The expression to compute the Deo index from 

Ipuxiey is 



Deo = 1-0488 (±0.0033) + 0.0051 (±0.0001) I Pux i ey 



(16) 
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Fig. 8. Comparison between measurements on the TNG subsample of the first CO bandhead using different index definitions. The 
panels show the transformation between the new index Deo and CO^ 8 , lFmgei> Ipuxiey and COkh (from top to bottom, and from left 
to right). The solid lines are least squares fits to the data, and correspond to the transformations given in Eqs. fT4l - [T7l 



with r 2 = 0.9629. 

Finally, we have also computed the conversion between the 
ratio COkh (not to be confused with CO^ g ; see Eq.[6]l and the 
new CO index 

D co = 2.1 119 (+0.0724)- 1.6205 (+0.1772) CO KH + 

+0.5521 (+0.1071) CO| H (17) 

with r 2 = 0.9867. This last transformation will be used in next 
section. 

4. Measurements of the CO absorption for the 
stellar library and error estimates 

4.1. Index measurements 

A detailed study of the continuum in spectra observed at CAHA 
compared with the the spectra published by KH86 have revealed 
some problems with the flux calibration of the CAHA data. The 
shape of the continuum in these spectra showed a spurious and 
non-reproducible high-frequency structure superimposed to the 
real continuum, which was affecting not only the shape of the 
continuum but also the final index measurements. During the re- 
duction of the data it was neither possible to identify nor correct 
this additional source of noise. 

In order to handle, at least in an empirical way, the spec- 
trophotometric calibration of the CAHA spectra, we re-observed 
a good subsample of the stellar library at the TNG. To guarantee 
that the TNG data were on the appropriate spectrophotometric 
system, the CO measurement of each star observed at the TNG 
was compared with the measurement of the most similar spec- 
trum (in r e ff and luminosity class) available in the KH86 library 
at the same spectral resolution. 



Since the bandpasses for the new index Dco encompass a 
wide range in wavelength, the strange behaviour of the contin- 
uum shape in the CAHA spectra has a large impact on the index 
measurements. Luckily, this is not such a big issue for the COkh 
ratio, since both continuum and absorption bandpasses are very 
close in this definition. For that reason, we decided to measure 
the COkh index, transforming afterward the results into the Dco 
index using Eq.[T7l( which provides a good conversion between 
both indices, as shown in the previous section). In more detail, 
the method to derive the final spectrophotometric calibration can 
be summarized as follows. 

First of all, the COkh measurements of the subsample of 
stars with solar metallicity re-observed at the TNG were com- 
pared with the corresponding star in the KH86 library, as ex- 
plained before. In Fig. [9] left panel, the results of this compar- 
ison are shown. A least squares fit to the one-to-one relation 
was computed, providing r 2 = 0.9768. Although some points 
in this figure appear relatively far from the 1 : 1 relation (consid- 
ering their error bars), it is important to highlight that we are not 
comparing the same stars, and that, in any case, the determina- 
tion coefficient r 2 is high enough to guarantee the quality of the 
transformation. 

Finally, we used the stars in common between the TNG 
and CAHA to empirically correct the measured indices on the 
CAHA spectra sampled at the spectral resolution of the TNG 
spectra. Fig. [9] right panel, presents the relation between the 
COkh ratio in common stars observed in both observatories. A 
least squares fit to a straight line gives r 2 = 0.9539. We have 
checked that a unique empirical correction is valid for all the 
observing runs at CAHA (within the statistical errors). For this 
study, we have used all the individual measurements for each 
star instead of averaging all the observations along the slit. After 
transforming the CAHA COkh measurements onto the correct 
spectrophotometric system, we applied Eq. [l7]to obtain the new 
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Fig. 9. Left panel: COkh index measured on the stars with solar metallicity observed at the TNG (CQ^ ) vs. the same measurements 
in the most similar stars from the KH86 library (COkh)- Right panel: COkh index measured on stars observed in common at TNG 
(CO£^ G ) and CAHA (CO^ HA ). Different symbols are used for distinct observing runs at CAHA, although no segregation of the 
data between runs is apparent. The dotted line shows the empirical correction for the CAHA measurements. In both panels, the 
dashed line shows the one-to-one relation. 



Deo measurements, which will be used later to derive the empir- 
ical fitting functions. 

4.2. Random errors and systematic effects 

4.2.1. Known sources of random errors 

We have considered several sources of random errors: photon 
statistic and read-out noise, errors due to the detector cosmetic, 
the combined effect of wavelength calibration and radial velocity 
uncertainties, and flux calibration uncertainties. 

(i) Photon statistics and read-out noise. With the aim of trac- 
ing the propagation of photon statistic and read-out noise, we 
followed a parallel reduction of data and error frames with the 
reduction package ReDjhE , which creates error frames at the 
beginning of the reduction procedure and translates into them, 
by following the law of combination of errors, all the manipula- 
tions performed over the data frames. In this way the most prob- 
lematic reduction steps (flat-fielding and distortion corrections, 
wavelength calibration, etc.) are taken into account and, finally, 
each data spectrum has its corresponding error spectrum which 
can be used to derive reliable photon errors in the index (cr photon ). 
A detailed description of the estimate of random errors in the 
measurement of classical line-st rength indices and the imp act of 
their interpretation are shown in lCardiel et a 1.(1998,2003). The 
new CO index is defined in this paper as a generic discontinuity 
and follows the expressions given in § l3.3l for the errors. 

(if) Errors due to the detector cosmetic. In the case of in- 
frared detectors, the errors due to the detector cosmetic are very 
important. We measured the r.m.s. {root-mean-squared) on two 
different regions of the spectra with no apparent absorption fea- 
tures OU2.2685-2.2790/jm and ^.2825-2.2890/im), in order 
to obtain an estimation of photon statistics and read-out noise er- 
rors together with the errors due to imperfections present in the 
images even after flat-fielding correction. In Fig.[l0]we compare 
the errors due to photon statistic and read-out noise errors (upper 
panel), derived from first principles (using the parallel reduction 
of data and error frames), with the errors from direct measure- 
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Fig. 10. Upper panel: Histogram with the photon statistic and 
read-out noise error in the Deo for the stars in the stellar li- 
brary. Lower panel: Similar histogram with the errors derived 
from measurements of the r.m.s. on two spectral regions with no 
apparent absorption features for the stars in the stellar library. 

ment of the r.m.s in the spectra (lower panel). As it can be seen, 
we are underestimating the random errors if we do not consider 
the errors due to the detector cosmetic. 

(Hi) Wavelength calibration and radial velocity errors. The 
combined effect of wavelength calibration and radial velocity 
correction is another random source of error. We considered as 
an upper limit for the wavelength calibration error the projection 
on the spectral direction of half of the slit width. These values 
were measured from the FWHM of the arc lines, providing er- 
rors of 45 km/s for the observations at CAHA and 85 km/s for 
the TNG images. These numbers translate into typical errors of 
cr[D co ] = 0.002 and cr[D co ] = 0.004, for CAHA and TNG 
respectively. 

Radial velocity for most of the stars in the stellar library wer e 
taken from the Hipparcos Input Catalogue (iTuron et al.lll992l) . 
which in the worst cases are given with mean probable errors of 
5 km/s. For the missing stars in this catalogue, we adopted radial 



velocities from the SIMBAD database, which presents typical 
errors of a few km/s. These errors are completely negligible in 
comparison with the wavelength calibration errors already men- 
tioned. 

(iv) Flux calibration errors. In order to check for possible 
random errors in the index measurements due to the flux cal- 
ibration method, we should have observed several spectropho- 
tometric reference spectra each observing night. Since this is 
a highly time-consuming approach, we did not observe them. 
Alternatively, we have decided to check the impact of this kind 
of error through the comparison of common stars between dif- 
ferent nights and runs and we will discuss it in § 14.2.31 

4.2.2. Additional sources of random errors 

Expected random errors for each star can be computed by adding 
quadratically the random errors derived from the known sources 
previously discussed, i.e., 



O 2 [CO] initial = 

= ^"^[COJphoton+cosmetic + °~ [CO] wavelength- 



(18) 



However, additional (but unknown) sources of random errors 
may still be present in the data. Using the multiple observations 
available for each particular star, we compared the standard de- 
viation of the different Deo measurements (cr[CO] rms ) with the 
initial random error (a 2 [Dcolimtiai) for that star. In the cases in 
which the standard deviation was significantly larger than the 
expected error (using a ;^ 2 -test of variances), a residual random 
error cr[Dco]residuai was derived as 



cr 2 [CO] residu ai = cr 2 [CO] rms - o 2 [CO] 



initials 



(19) 



and quadratically added to the initial random errors to obtain the 
final random errors, i.e., 



cr 2 [CO] fina i = cr 2 [CO] init iai + o 2 [CO] 



residual ■ 



(20) 



It is worth noting that this additional error was only nec- 
essary for a few stars, since most of them presented consistent 
errors between different measurements. 



4.2.3. Systematic effects 

Due to the large number of runs needed to complete the whole li- 
brary, systematic errors can arise because of possible differences 
between the spectrophotometric systems of each run at both the 
CAHA and TNG telescopes. To guarantee that the whole dataset 
is on the same system, we observed stars in common between 
different runs at each telescope. We compared the index mea- 
surements of these stars by applying a f-test to study whether 
the differences of the index measurements of the common stars 
were statistically larger than zero, finding no systematic effects 
between different nights within a given observing run, and be- 
tween different runs. For that reason we have considered that 
an additional error exclusively due to flux calibration was not 
required (i.e., the actual flux calibration errors are within the al- 
ready computed random errors). 

It is important to highlight that there is not guarantee that 
our data are on a true spectrophotometric system, so we encour- 
age readers interested in using our results to include in their ob- 
serving plan stars in common with the stellar library to ensure a 
proper comparison. 



4.3. Additional measurements of the CO absorption: globular 
cluster stars 

To improve the stellar parameter coverage of our stellar library, 
additional stars were included for the co mputation of e mpiri- 
cal fitting functions for the D eo (see § [6j. Fro gel et all {2001) 
and IStephens & Fro gel (2004) presented a sample of globular 
cluster giant stars (R ~ 1500, <x ~ 85 km/s), characterized by 
their low metallicity, with mea surements of the CO absorption 
at 2.3 //m using the definition of Froge feTaTI (1200 ll) . Since there 
is not dependence of iFrogei at the spectral resolution of these 
data (see Fig [5]), we transformed their CO measurements to the 
new index by applying Eq. [15] The stellar atmospheric param- 
eters of these stars were determined from J and K photometry, 
as it is explained i n § 15.21 Finally, w e considered 80 stars from 
IFrogei etal](l200lh and 14 stars from lStephens & Frogell(l2004 . 
which, together with the stellar library presented in this work, 
will be used to parametrize the behaviour of the CO index as a 
function of the stellar atmospheric parameters. 



5. Stellar atmospheric parameters 

In this section, we present the atmospheric parameters for the 
stars considered in the computation of the empirical fitting func- 
tions for the Deo index presented in § [6j 



5.1. Atmospheric parameters for the observed stellar library 
sample 

As we mentioned before, the stellar library presented 
in this work is a subsample of the MILES library 
( Sanchez-Blaz quez et alj I20 06). Following the previous 
work by ICenarro et alJ (l2001bl) . a reliable and highly homoge- 
neous data set of atmospheric stellar parameters fo r the stars in 
MILES library was derived by Cenarro et al. I d2007l) as the result 
of a previous, extensive compilation from the literature and a 
carefully calibration and correction by bootstrapping of the data 
on to a reference system . In short, the method can be itemized in 
the fo llowing steps (see Cen arro et al.ll2001bl and Cenarr o~et al.l 
2007 for a more detailed explanation of the working procedure): 
(i) selection of a high-quality, standard reference of atmospheric 
parameters, (ii) bibliographic compilation of atmospheric pa- 
rameters of the stars in the library, (iii) calibration and correction 
of systematic differences between the different sources and the 
standard reference system, and (iv) determination of averaged, 
final atmospheric parameters for the library stars from all those 
references corrected on to the reference system. Together with 
the atmospheric parameters, an estimation of the uncertainties 
in their determination is given for each star. In that way, stellar 
atmospheric parameters (and their uncertainties) are perfectly 
well known for the stars in the new stellar library. 



5.2. Atmospheric parameters for the additional globular 
cluster stars 

For the subsample of 94 red giant branch (RGB) stars from 
Galactic globular clusters, effective temperatures (T e s) and sur- 
face gravities (logg) were derived from J and K photometry, 
follow ing a tec hnique similar to that explained in ICenarro et al.l 
d2007l) (see also lCenarro et alj|2001bl) . Absolute, re ddening cor- 
rected pho tometry for all clus t ers w as taken from Frog el etafl 
(I2001I) and lStephens & Frogel (I20041) . 

The surface gravity for each cluster star was estimated by 
matching its location in a Mk-( J - K) diagram to the appropri- 
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Fig. 11. Pseudo-HR diagrams for the cluster stars together with adequate isochrones dGirardi et al.ll2000tlLejeune et al.lll997l. fl 998) 
for each individual cluster. Open circles are used for individual stars in the clusters. Solid and dashed lines illustrate isochrones of 
12.59 Gyr and two metallicity values enclosing the one of the cluster (as shown in the labels). Adopted metallicities for the clusters 
are labeled in the figure. In all cases, the solid line is employed to indicate the isochrone whose metallicity is closer to that of the 
cluster. Surface gravity for each star was estimated by comparing to the predicted Mk as explained in § 15.21 Final atmospheric 
parameters for each cluster star, and their corresponding errors, are given in Table I A. 21 



ate isochrones from iGirardi et al.l d2000t) (see Fig. [TTT i. whose 
colors and magnitudes were previously transformed to the ob- 
servational plane using the empirical colour -te mperature rela- 
tions for gi ant stars from lAlonso et al.l (1 19991) and lLeieune et al.l 
(1 19971. 1 19981) . with the latter being for gi ants with r eff < 3500 K 
(see more details in Vazdekis et al. 20031). 

If availab le, the metallici t y for each cluster was taken from 
the work by iRutledge et all (119971) (hereafter RHS97), which 



provi des metallicity es t imate s for Galactic globular cluster on 
the ICarretta & Grattonl (1 19971) (hereafter CG97) scale based on 
the Call triplet strengths of their RGB stars. This was the case 
for NGC0104, NGC0288, NGC0362, NGC5927, NGC6553, 
NGC6624, NGC6712, and M69. F or M7 1 , however, we kept the 
value in the CG97 scale inferred by lCenarro et al ] (120021) accord- 
ing to the CaT indices of their stars. If not av ailable in RHS97, 
metallicity values in the (IZinn & West! 11984 hereafter ZW84) 
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Fig. 12. Behaviour of the first CO bandhead with T e ff, logg and [Fe/H] (from left to right). In each case, the displayed stars have 
been selected to exhibit a range in a given atmospheric parameter (as shown in the labels) with similar values for the other two 
parameters. 



scale were transformed into the CG97 scale using Equation 7 
in CG97. This was the case for the rest of our GC sample. For 
NGC6388, NGC6440, Lillerl, and Terzan2, ZW84 metallicity 
values were taken from Table 6 in that paper. For NGC6528, 
the value given in the 2003 revised version of the catal og of 
parameters for Milky Way globular clusters dHarrisl ll996h was 
employed. It is important to note that there is a double reason 
for using the CG97 metallicity scale in this work. First, since, 
as compared to the IZinn & Westl (1 19841) metallicity scale, the 
agreement between the locus of the cluster RGB stars and the 
corresponding isochrones in the Mk-(J - K) plane is much bet- 
ter when using CG97 values, particularly in the high metallicity 
regime. Also, because the metallicit ies of the rest of cluster stars 
in this paper (as taken from either ICenarro et al" I l2001bl 120071) 
are already on the CG97 scale, thus guaranteeing full consis- 
tency among the metallicities of the whole star sample, and min- 
imizing systematics in the computation of the fitting functions. 

Assuming that all the galactic globular cluster under study 
are similarly old, and taking into account the above metal- 
licities, for all the stars in each cluster we used two of the 
Girardi's isochrones of 12.6 Gyr with metallicities enclosing the 
corresponding value of the cluster (solid and dashed lines in 
Fi g . fTTb . For each isochrone, a \ogg value for each star was es- 
timated by interpolating in Mk- Final logg values were com- 
puted as weighted means of the single values derived from the 
two different metallicity isochrones, with weights accounting for 
the differences between the adopted cluster metallicity and the 
isochrone values. 

Effective temperatures for all the cluster giants were com- 
puted on th e basi s of the (J - K)-T e s relati ons given in 
lAlonso et alJ d 19991) (for T eff > 3500 K) and iLeieune et all 
(119971 119981) (for T eff < 3500 K). Interestingly, since Alonso's 



relation for (J - K) does not depend on either metallicity and 
log g, T e g can be directly determined from (J - K). This indeed 
helps to minimize the sources of uncertainties in the final tem- 
peratures of most globular cluster giants. As a matter of fact, 
only a few temperatures for very cold stars were computed us- 
ing Lejeune's relations. 

Errors for the derived atmospheric parameters of the cluster 
stars were estimated from input uncertainties in [Fe/H] and in the 
photometric data employed in this technique . For t hose [Fe/H] 
values taken from RHS97 and lCenarro et alJ d2002h . the uncer- 
tainties given therein were assumed. For the rest of the clus- 
ters, metallicity errors were computed from those given in ZW84 
through stan dard error propag ation in E quation 7 of CG97 . Since 
most stars in lFrogel etail d200ll) and lStephens & Fro gel (2004) 
were selected from the bright ends of the globular cluster lu- 
minosity functions, photometric observational errors in J and 
K bands turned o u t to be very small (« 0.06 mag; see e.g. 
Frog el et alJ 119951; iKuchinski et ail 119951; IKuchinski & Frogell 
1995) as compared to typical uncertainties in the assumed dis- 
tance moduli and reddening corrections. Errors in Mk, and 
( J - K) are therefore strongly dominated by the above effects. 
For this reason, typical errors of 0.20 mag and 0.06 mag for Mk 
and ( J - K), respectively, have been considered for all the clus- 
ter stars. 

Final values of T e g, \ogg, [Fe/H], and their corresponding 
errors for the whole sample of cluster stars are presented in 
Table|A2] 

6. Empirical calibration of the new CO index 

In this section we parametrize the behaviour of the new CO in- 
dex in terms of the stellar atmospheric parameters (T e ff , log g and 
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Fig. 13. D co index as a function of (= 5040/r eff ) for the 
stars of the new l ibrar y plus the globular cluste r stars from 
IFrogel et all (l200ll) and IStephens & Frogell (12004 used in this 
study. Symbols mean metallicity, as explained in Fig.[TJ Relative 
symbol sizes indicate gravity ranges (small symbols for dwarfs, 
increasing symbol size for decreasing gravity). 



[Fe/H]. For that purpose, we use the Dco index measurements of 
the stars of the new stellar library and o n th e sample of globu- 
lar clu ster stars from lFrogel et al.l (1200 ll) and Stephens & Frogel 
(2004) (see § O- This paper shows a detailed quantitative 
metallicity dependence of the CO feature at 2.3 pm. 

6.1. Behaviour of the CO index as a function of stellar 
atmospheric parameters 

The dependence of the strong CO absorption bands as a function 
of the effective temperature and surface gr avity is well known 
from t he first studies in the infrared (KH86, Johnson & Mendez 
1970t IFrogel et all [l978t lLancon & Rocca-Volmerangelj 199 ' 
Dovon et al.ll99 4: Wallace & Hinkle 1997; Ramirez et al.ll99 



Forster Schreiberl 120001: lLancon & Wood! l2000t llvanov et all 
2004; ISilva et al.ll2008l) . The CO absorptions turn deeper as the 
stars are cooler, while hot stars show no trace of CO features. 
On the other hand, as surface gravity decreases, the CO ab- 
sorptions become prominent, i.e., supergiant stars present more 
important CO absorption than dwarfs. All the previous works 
in the K band highlighted both the tight correlation of their 
defined CO indices with temperature (spectral type or J - K 
color in the first papers), and the dependence of the CO absorp- 
tion with surface gravity. In addition, a few theo retical studies 
dMcWilliam & Lambed 1 19841 lOridia et al.lll993l) indicate that 
these spectroscopic f eatures should be metallicity dependent. 
iTerndrupet all (1199 ll) showed that Baade's window stars have 
deep 2.3 pm bands, what they interpreted as these s tars prob- 
ably h aving a high metallicity. Model predictions by iMarastonl 
(2005) also showed the dependence of the CO absorptions in 
the K band with metallicity. Observations of composite s tellar 
population (e.g. lOriglia & Olival 120001: iRiffel et allfeoOTl) give 
support to this dependence. Unfortunately, the lack of empirical 
stellar libraries with an appropriate coverage in metallicity had 
prevented, until now, a detailed quantitative description of this 
dependence. Fig. IT2l shows, graphically, how the first CO band- 
head at 2.3 pm varies with the stellar atmospheric parameters 
(T e s, logg and [Fe/H]). 



Several authors have partially parametrized the described be- 
haviou r of the first CO bandhead. lLancon & Rocca-Volmer ange 
(1992) computed a relati on of their CO ind ex with the color 
temperature of giant stars. iDoyon et al.l d 19941) parametrized the 
behaviour of their CO index with the effective temperature for 
supergiant, gi a nt and dwarf stars, separately. More recently, 
Ram irez et al.l (1 19971) used their CO index to obtain effective 
temperatur es for gian t s. The y applied the same method to dwarf 
stars from lAli et alj ( 1 1995b . However, there is no systematic 
study of the dependence of the CO absorption as a function of 
the atmospheric stellar parameters due to the poor stellar param- 
eter coverage of previous libraries, especially in metallicity. In 
the next subsection we compute an empirical calibration for the 
Dco index measured for the stars of our stellar library which ac- 
counts for the previously described qualitative behaviour of the 
CO absorption. 

6.2. Fitting functions 

To reproduce the behaviour exhibited by a given feature, em- 
pirical calibrations of the corresponding line-strength index as 
a function of the stellar atmospheric parameters were derived. 
These calibrations, called fitting functions, are just a mathemat- 
ical description of the observed behaviour and there is not phys- 
ical justification for the explanation of each of the significant 
terms in such fitting functions. We use 6 = 5040/r e ff as temper- 
ature indicator, together with log g and [Fe/H] as p arameters for 



gravity and metallicity. Foll o wing previous work s (Gorgas et al 
' IWorfheyetaDll994t iGorgas et al.l 119991: ICenarro et al 



1993 



2002), two possible functional forms for the computation of the 



fitting functions are typically considered, in particular 
/ (6, log g, [Fe/H]) = p (0, log g, [Fe/H]) 



(21) 



and 



1(8, logg, [Fe/H]) = constant + exp[p (0, \ogg, [Fe/H]) , (22) 

where / refers to a given index, and p (6, log g, [Fe/H]) is a poly- 
nomial with terms up to the third order, including all possible 
cross-terms among the three parameters, i.e. 

p (0, log g, [Fe/H]) = J] c iJtk ff (logg) J [Fe/U] k , (23) 

()</+ ;'+it<3 

with < i + j + k < 3 and < i, j, k. 

The optimum fitting function is the one which minimizes the 
residuals of the fits, i.e., when the differences between the mea- 
sured index (7 bs) and the index predicted by the fitting function 
(7 pre d) are the lowest. 

In general, not all the terms in Eq. [23] are necessary for re- 
producing the dependences of an index on the stellar parameters. 
To obtain the significant terms fo r the best fitting func tions, we 
follow two different approaches (Cen arro et al.l 12002). Both of 
them consist of an iterative and systematic procedure based on 
the computation of a general fit together with the analysis of the 
residual variance of the fit and estimated errors of each fitted co- 
efficient, for a given set of stars with well known atmospheric pa- 
rameters, index measurement and error. The significance of each 
term considered in the fit is calculated by means of a f-test (that 
is, using the error in that coefficient to check whether it is sig- 
nificantly different from zero). Typically, we consider a term as 
significant for a threshold value of the significance level parame- 
ter a = 0.10. The first method consists of computing the fit with 
all the terms in Eq. [23] and eliminating, during each iteration, 



the least significant term. A parallel method consists of comput- 
ing an initial one-parameter fit and incorporating new terms only 
when they are significant. Together with this criterion, we study 
the residuals derived for each set of terms. The final fit will be 
the one which minimizes the sum of the residuals while having 
all the employed c, ^ coefficients statistically significant, that 
does not produce systematic deviations in the residuals for stars 
of different types, e.g., stars of a given metallicity range, and that 
reproduces the observed behaviour of the index. 



Besides the global behaviour described for giant stars, two 
different trends are found for this type of stars around 9 = 
1.3 - 1.4 in Fig. Qj] After a carefully study of these stars with a 
higher CO index, we found that they are stars in the asymptotic 
giant branch (AGB). For that reason, we decided to compute an 
independent fit for these stars in the range 9 — 1 .01 - 1 .56, shown 
Fig-ERcoefficient in Tabled Since there are no AGB stars for 
9 > 1 .56, we simply extrapolate the constant value of the CO 
index at 9 = 1.56. 



6.2.1. The general fitting procedure 

As a result of the large stellar parameter coverage of the library, 
there is not a single fitting function able to reproduce the whole 
behaviour of the Deo index. For that reason, we have divided 
the stellar parameter space into several ranges (boxes) where lo- 
cal fitting functions have been computed following the methods 
explained in the previous section. The final fitting function for 
the whole parameter space is then constructed by considering 
the derived local fitting functions. In some ranges, an overlap- 
ping region in a generic parameter z has been allowed between 
two different boxes. If I\ (x,y,z) and 1% (x,y,z) are two local fit- 
ting functions defined respectively in the intervals {z\,\,Z\,2) and 
(Z2,\,Z2,2) with Zz,i < Zi,2> the final predicted index I(x,y,z) in 
the overlapping region is obtained by interpolating the index 
value in both boxes, i.e., 



I(x,y,z) = wl\ (x,y,z) + (1 - w)h (x,y,z) ■ 



(24) 



where the weight w is modulated by the distance to the overlap- 
ping limits as 



w = 



Z2,l 



Zl,2 - £2,1 



or 



w = cos 



/ Z-Z2A 



(25) 



(26) 



_2 \Zi,2 - Z2.1 1 

for a cosine-weighted mean, where zi,\ <z<Z\,i- 

6.2.2. Fitting functions for the D C o index 

Fig. Q~3] shows the behaviour of the Deo index versus 9 for the 
stars in the new library. There is a clear dichotomy in the be- 
haviour of stars depending on their gravity. For that reason, 
we have divided the stellar atmospheric parameter space in two 
main groups: dwarf stars (logg > 3.5 dex) and giant and super- 
giant stars (logg < 3.5 dex). As we explained in § 16.11 there is 
also a strong dependence of the CO absorption with the effective 
temperature. That is why we have subdivided each gravity group 
into different temperature ranges. Independently of their gravity, 
stars with high effective temperature exhibit no traces of CO ab- 
sorptions and their index value tends to a constant (Deo - 1.05). 
On the other hand, due to the lack of very cold stars in both grav- 
ity regimes, we have computed a constant value of the index for 
cold dwarf and giant stars. In short, we have considered three 
temperature ranges for dwarf stars, while we have considered 
four different ranges for the giants (see Tabled. The boundaries 
of these ranges are shown in Fig. Q] In some cases, an overlap- 
ping region has been considered. After some experimentation in 
order to obtain the smoothest fit, we have computed the final in- 
dex in the intersection region considering a different weight w 
depending on the case (see Tabled. 



Up to now, we have described the general method to com- 
pute fitting functions. However, we have used a slightly different 
method for cold and warm giant stars. First of all, we have ob- 
tained constant fits for hot (9 < 0.90) and cold stars (8 > 1.55) 
as explained before. After this, we have generated a set of fake 
stars with random gravity and metallicity, and effective tempera- 
ture of 9 = 5040/r e ff = 1.60. Their index values were evaluated 
from the fitting function computed for cold giant stars. For the 
associated error of the indices of these fake stars we assumed the 
mean error of the cold stars. We calculated the local fitting func- 
tion of cool giants, including these new stars. In this way, we 
were forcing the fitting function to pass through these stars, i.e., 
approaching the constant value for cold giants. In a similar way, 
for the computation of the fitting functions of warm giants, we 
also created fake stars at 9 — 1.10 and 9 = 1.13, and we assigned 
their index value according to the fitting function for cool stars. 
A new set of random stars at 9 — 0.92 was introduced with the 
constant index value computed from hot stars. Finally, we used 
these fake stars together with the real stars to obtain the final fit 
for the warm giants. 

We derived the local fitting functions computing a weighted 
least squares fit to the stars within each parameter range, with 
weights according to the inverse of the squared uncertainties 
of the indices for each individual star. Here we considered the 
nominal errors in the index measurements, using the uncer- 
tainties derived in § 14.21 for the librar y stars, and t he ind ivid- 
ual errors for each clust er quoted by iFrogel et al.l (1200 ll) and 
Stephens & Frogel (2004). After the analysis of the residuals, 
it was necessary to repeat the fitting procedure with the inclu- 
sion of additional random uncertainties in some of the stars (see 
next subsection). The final fitting functions for each stellar pa- 
rameter range are presented in Table [6] This table includes the 
functional forms of the fits, the significant coefficients and their 
corresponding formal errors, the typical index error for the N 
stars employed in each interval (crj yp = N/ £j=i err 1 ), the un- 
biased residual variance of the fit (cr 2 ) and the determination 
coefficient (r 2 ). Also the intersection regions are indicated and 
the procedure used for computing the index. We plot in Fig. [14] 
the fitting functions for giant and dwarf stars as a function of 9 
for several metallicities, as well as the simple fit for AGB stars. 
Fig. [T5l shows the final residuals (ADco = Dcoobs - Dcopied) of 
the new CO index for the different groups of stars used for the 
computation of the fitting functions (stars observed at CAHA, 
stars observed at the TNG and cluster stars from previous works) 
as a function of the atmospheric parameters. As we expected, the 
residuals for each group of stars do not exhibit any systematic 
trend. 

Users interested in implementing these fitting functions 
into their population synthesis codes can make use of 
the FORTRAN subroutine available at the URL address 
http : //www . ucm . es/inf o/Astrof /ellipt/CO . html 




©=(5040/T eff ) 

Fig. 14. Empirical fitting functions for the Deo index for giants, including the fit for AGB stars (top panel) and dwarfs (bottom 
panel), computed as explained in the text. Symbol types and size mean metallicity and gravity, the same than in previous figures. 
The different lines represent the metallicities [Fe/H]= +0.5, +0.0, -0.5, -1.0, -1.5, and - 2.0, from top to bottom. The reduced 
scatter for 9 < 0.9 is due to the absence of CO absorptions for hot stars, as explained in the text. 



Table 6. Coefficients and statistical data for the local empirical 
fitting functions for the Deo index in each range of stellar pa- 
rameters. 



Hot dwarfs 


0.38 < 9 < 0.90 


3.50 < logg < 5.50 


exponential fit 


N = 28 stars 


cr typ = 0.00913 




0.0499 ±0.0011 


o- std = 0.00557 


Intersection 


Cosine-weighted mean 


0.80 < 9 < 0.90 


Cool dwarfs 


0.80 <0< 1.50 


3.50 < logg < 5.50 


exponential fit 


N = 39 stars 


cr typ = 0.00748 


Co 


-0.0292 ± 0.0330 


cr std = 0.01254 


9 


0.1006 ± 0.0329 


r 2 = 0.765 


[Fe/Hl 


0.0174 ± 0.0067 




Intersection 


Cosine-weighted mean 


1.45 < 6 < 1.50 


Cold dwarfs 


1.33 < 9 < 1.80 


3.50 < logg < 5.50 


exponential fit 


N = 7 stars 


cr typ = 0.00960 


Co 


0.1025 + 0.0046 


cr std = 0.01195 


Hot giants 


0.42 < 9 < 0.90 


-0.40 < logg < 3.50 


exponential fit 


N = 15 stars 


cr typ = 0.01954 


Co 


0.0459 ± 0.0010 


o- std = 0.00398 


Intersection 


Cosine-weighted mean 


0.90 < 9 < 0.93 


Warm giants 


0.90 < 9 < 1.131 


0.13 logg < 3.5 


exponential fit 


N = 63 stars 


cr typ = 0.00764 


CO 


-0.3073 ± 0.0046 


cr std = 0.00428 


e 


0.3876 ± 0.0043 


r- = 0.982 


[Fe/H] 


-0.1016 ±0.0042 




0[Fe/H] 


0.1072 ±0.0039 




[Fe/H] 2 


-0.0023 ± 0.0005 




Intersection 


Distance-weighed mean 


1.09 < e < i.io 


Cool giants 


1.10 <9< 1.60 


-0.34 < logg < 3.41 


exponential fit 


N = 167 stars 


cr, yp = 0.01062 


Co 


-0.5224 ± 0.0970 


cr std = 0.00890 


6 


0.8257 ±0.1417 


r 2 = 0.958 


[Fe/H] 


0.0674 ± 0.0101 




0[Fe/H] 


-0.0444 ± 0.0065 




e 2 


-0.2200 ± 0.0509 




[Fe/H] 2 


-0.0023 ± 0.0014 




No intersection 


Cold giants 


1.55 < 9< 2.03 


-0.07 < logg < 3.50 


exponential fit 


N = 7 stars 


cr typ = 0.02156 


Co 


0.2397 ± 0.0107 


o- std = 0.02698 


AGB stars 


1.01 <9< 1.56 


-0.11 <logg< 1.56 


exponential fit 


N = 18 stars 


o- lyp = 0.00612 


Co 


-0.8893 ± 0.2198 


cr std = 0.00892 


6 


1.4950 ±0.3610 


r 2 = 0.985 


6- 


-0.4816 ±0.1461 





6.2.3. Residuals and error analysis 

To explore in more detail the reliability of the fitting functions, 
we computed the unbiased residual standard deviation from the 
fits, o- s id = 0.0093, and the typical error in the measured index, 
o"typ = 0.0025, for the stars employed in the computation of the 
global fitting functions derived above. 

After the initial fit, we discovered that <x std was larger than 
what should be expected uniquely from index uncertainties (see 
also the values of cr st d and <x typ for different subgroups of stars 
in Table [7]). Furthermore, a F test of comparison of variances 
showed that <x std was statistically larger than cr typ , leading to the 
introduction on an additional residual error (<x;? es = cr^ [d - of yp ). 
Such residual error could be due to uncertainties in the input 
stellar atmospheric parameters, not included in the error bud- 
get so far. In order to check this, we have computed how errors 
in the input stellar parameters translate into uncertainties in the 
predicted Deo- This depends both on the local functional form 
of the fitting function and on the atmospheric parameters range 



Table 9. Absolute errors in the fitting functions predictions for 
different values of the atmospheric parameters. Input logg val- 
ues varying with effective temperature for dwarfs (V), giants 
(III) and supergiants (I) have been taken from lLangl (1 19911) . 



Teg 


[Fe/H] 


V 


o-[D co ] 
III 


I 


6000 




0.002 


0.001 


0.001 


5500 


+0.5 


0.010 


0.002 


0.002 


5500 


+0.0 


0.007 


0.001 


0.001 


5500 


-1.0 


0.006 


0.001 


0.001 


5500 


-2.0 


0.011 


0.002 


0.002 


5000 


+0.5 


0.010 


0.001 


0.001 


5000 


+0.0 


0.007 


0.001 


0.001 


5000 


-1.0 


0.005 


0.001 


0.001 


5000 


-2.0 


0.011 


0.001 


0.001 


4500 


+0.5 


0.011 


0.004 


0.004 


4500 


+0.0 


0.008 


0.003 


0.003 


4500 


-1.0 


0.007 


0.003 


0.003 


4500 


-2.0 


0.011 


0.006 


0.006 


4000 


+0.5 


0.014 


0.003 


0.003 


4000 


+0.0 


0.012 


0.002 


0.002 


4000 


-1.0 


0.011 


0.003 


0.003 


4000 


-2.0 


0.014 


0.005 


0.005 


3500 


+0.5 


0.018 


0.003 


0.007 


3500 


+0.0 


0.018 


0.003 


0.007 


3500 


-1.0 


0.018 


0.003 


0.007 


3500 


-2.0 


0.018 


0.004 


0.007 


3000 




0.005 


0.014 


0.014 



(e.g. hot stars have T e f[ uncertainties larger than cooler stars). 
For each star of the sample, we have derived the errors cor- 
responding to the uncertainties in effective temperature (cr^), 
gravity (o-\ ogg ) and metallicity (cr [Fe/H])- The effect of the un- 
certainties in the three stellar parameters were finally computed 
as crJL = cr\ + cr? + cr}^ Table 17] presents, summarized 
for the different stellar groups, all the error estimations already 
discussed. Finally, in the cases where the residual errors (cr res ) 
were not explained by the uncertainties in the stellar parameters 
(<x pal ), an extra residual error was added to the original index er- 
ror. At the end, no iterations were needed for the globular cluster 
stars and AGB stars, one iteration was required for the giants 
and two iterations for dwarfs (in any case, the additional error 
for dwarf stars is lower than the necessary for giants). The un- 
certainties of the final Deo fitting functions are listed in Table[8] 
Note that although in this table cr std is still larger than cr typ for 
the four initial subgroups of stars analyzed in Table|7j both stan- 
dard deviations are statistically comparable using the F test of 
variances previously mentioned. 

Finally, as the purpose of this paper is to predict reliable in- 
dex values for any combination of input atmospheric parameters, 
we have also computed the random errors in such predictions. 
These uncertainties are given in Table |9]for some representative 
values of input parameters. 

7. Summary 

The aim of this work was to obtain an accurate empirical calibra- 
tion of the behaviour of the CO feature at 2.3 fim for individual 



Table 7. Uncertainties of the initial Deo fitting functions for different groups of stars and mean Deo errors due to uncertainties in the 
input atmospheric parameters. N: number of stars. <x std : unbiased residual standard deviation for the fit. cr typ : typical observational 
Deo error for the subsample of stars. cr les : residual error. <jt cS and o~ [Fe/H] • mean Deo errors due to uncertainties in the input T e g 
and [Fe/H] (no error due to uncertainties in the input log g is computed since the computed fitting functions do not require gravity 
terms). cr par : total error due to atmospheric parameters (quadratic sum of the previous errors). <x std is not explained by the cr typ and 
cr par values for giants and dwarfs, and a residual error has been considered for these group of stars 





N 


0"std 


0"typ 


0"res 




f"[Fe/H] 




Field dwarfs 


54 


0.0059 


0.0017 


0.0056 


0.0010 


0.0009 


0.0014 


Field giants 


147 


0.0108 


0.0023 


0.0105 


0.0062 


0.0019 


0.0065 


Globular cluster stars 


85 


0.0276 


0.0236 




0.0139 


0.0010 


0.0140 


AGB stars 


19 


0.0115 


0.0061 




0.0097 


0.0000 


0.0097 


All 


305 


0.0093 


0.0025 




0.0075 


0.0014 


0.0078 



Table 8. Uncertainties of the final Deo fitting functions for different groups of stars and same computed errors than in Table [7] 





N 


0"std 


°~typ 


°> Efr 


O" [Fe/H] 




Field dwarfs 


54 


0.0086 


0.0078 


0.0010 


0.0013 


0.0017 


Field giants 


147 


0.0123 


0.0113 


0.0062 


0.0014 


0.0064 


Globular cluster stars 


85 


0.0252 


0.0236 


0.0135 


0.0009 


0.0135 


AGB stars 


19 


0.0115 


0.0061 


0.0097 


0.0000 


0.0097 


All 


305 


0.0130 


0.0114 


0.0074 


0.0012 


0.0076 
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Fig. 15. Residuals of the derived fitting functions (observed minus derived) against the th ree input stel l ar atm osp heric parameters for 
the sta rs observed at CAHA (upper panels) and TNG (middle panels). Cluster data from Frogel et al. (200lJ) and lStep hens & Frogel 
(2004) are presented in lower panels. Codes and relative sizes of the star symbols are explained in Fig. [TJ and [13] 



stars, with the purpose of making it possible to obtain reliable 
predictions for the CO strength for stellar populations in unre- 
solved systems with a wide range of ages and metallicities. The 
main results of this work can be summarized as follows: 

1. We present a new stellar library in the spectral re- 
gion around the first CO bandhead at 2.3 fim. It con- 
sists of 220 stars with stellar atmospheric parameters in 



the range 2485 < r eff < 13404 K, -0.34 < logg < 5.30 dex, 
-2.63 < [Fe/H] < 0.98 dex. 
2. We define a new line-strength index for the first CO band- 
head at 2.3 fim, Deo, l ess sensitive to spectral resolution, 
wavelength calibration, signal-to-noise ratio and flux calibra- 
tion than previous definitions. 



3. We compute empirical fitting functions for the Deo to 
parametrize the behaviour of the CO feature as a function of 
the stellar atmospheric parameters. In this work we quantify, 
for the first time, the metallicity dependence. 

We expect that the work presented in this paper will help 
researchers to start exploiting in depth the so far poorly- 
explored and poorly-understood near-IR spectral region centered 
at 2.3 yum, since, as we have shown, the strong CO bandhead can 
be employed to extract useful physical information of composite 
stellar populations. 
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Appendix A: Stellar library and cluster stars 

In Table lA. fl we list the stars of the stellar library with their spec- 
tral type, K magnitude, effective temperature and the associated 
uncertainty (T e ff and <x[T e jf]), surface gravity and its uncertainty 
(logg and <r[logg]), metallicity and its uncertainty ([Fe/H] and 
cr[Fe/H]), number of measurements (N), Deo index and their er- 
ror cr[Dco]- The stars observed at the TNG are labeled with a | 

in the index measurement. 

In Table lA. 2l we list the cluster stars from Frogel et al. (2001) 
and Stephens & Frogel (2004) employed in the computation of 
the empirical fitting functions. We present the derived atmo- 
spheric parameters and their associated uncertainties, and the 
Deo index and error for each star. 



In both tables, AGB stars are labeled with a ★ in the name of 
the star. 



Table A.l. Stellar library used in the fitting function procedure. Stars marked with a f symbol were observed at the TNG. AGB stars are marked 
with a ★ symbol. 
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0.50 


-0.93 


0.30 


3 


1.057 


0.009 


HD078712 


M6SI 


-1.87 


3202 


60.9 


0.00 


0.18 


-0.11 


0.09 


4 


1.216* 


0.010 


HD078732 


G8II 


+3.20 


4900 


117.6 


2.00 


0.21 


+0.24 


0.10 


8 


1.108 


0.011 


HD079211 


MOV 


+4.14 


3710 


60.9 


4.71 


0.18 


-0.40 


0.10 


3 


1.093 


0.010 


HD079452 


G6III 


+3.86 


4829 


60.9 


2.35 


0.18 


-0.84 


0.09 


4 


1.086* 


0.012 


HD081192 


G7III 


+4.11 


4705 


75.0 


2.50 


0.40 


-0.62 


0.15 


4 


1.101* 


0.011 


HD082074 


G6IV 


+4.15 


5055 


117.6 


3.30 


0.21 


-0.48 


0.10 


2 


1.088 


0.011 


HD082885 


G8IV-V 


+3.69 


5488 


60.9 


4.61 


0.18 


+0.00 


0.09 


2 


1.061 


0.009 


HD083425 


K3III 


+ 1.58 


4120 


60.9 


2.00 


0.18 


-0.35 


0.09 


4 


1.170* 


0.011 


HD083618 


K3III 


+0.87 


4231 


60.9 


1.74 


0.18 


-0.08 


0.09 


4 


1.161* 


0.011 


HD083632 


K2III 


+4.72 


4214 


60.9 


1.00 


0.21 


-1.39 


0.09 


4 


1.160* 


0.011 


HD085235 


A3IV 


+4.37 


11200 


117.6 


3.55 


0.21 


-0.40 


0.10 


3 


1.049 


0.009 


HD085503 


KOIII 


+ 1.36 


4472 


75.0 


2.33 


0.40 


+0.23 


0.15 


3 


1.136 


0.011 


HD085773 


G:w?I 


+7.95 


4463 


60.9 


0.98 


0.18 


-2.17 


0.09 


1 


1.056 


0.014 


HD087737 


AOIb 


+3.29 


9625 


60.9 


1.98 


0.21 


-0.04 


0.10 


4 


1.053* 


0.010 


HD087822 


F4V 


+5.13 


6590 


60.9 


4.15 


0.18 


+0.14 


0.09 


4 


1.049* 


0.009 


HD089484 

X XXV V/O y i U i 


KHIIb 


-0.81 


4470 


60.9 


2.35 


0.18 


-0.49 


0.09 


3 


1.119 


0.011 


HD089822R 


AOsn^TTT 


+3.39 


5538 


117.6 


2.44 


0.21 


+0.51 


0.10 


2 


1.050 


0.010 


HD092523 


K3III 


+ 1.55 


4090 


60.9 


1.96 


0.18 


-0.38 


0.09 


2 


1.163 


0.011 


HD093487 


F8 


+6.76 


5250 


60.9 


1.80 


0.18 


-1.05 


0.09 


2 


1.068 


0.012 


HD095578 


MOIII 


+0.80 


3700 


60.9 


1.40 


0.18 


-0.23 


0.09 


2 


1.206 


0.011 


HD095735 


M2V 


+3.25 


3551 


60.9 


4.90 


0.21 


-0.20 


0.09 


8 


1.108* 


0.009 


Hn09fS^fiO 




+2 If, 




1 1 7 fi 


so 


21 


—0 ss 
w — >o 


1 

U. 1U 


A 


1 244"'" 




HD097907 


K3III 


+2.43 


4351 


60.9 


2.07 


0.18 


-0.10 


0.09 


2 


1.162 


0.011 


HD099648 


G8II-III 


+2.83 


4850 


117.6 


1.90 


0.21 


+0.36 


0.10 


3 


1.097 


0.011 



Table A.l. continued. 



Name 


Spectral 


K 


/ ofF 

* err 


Tp.fi\ 
u L 1 err j 


' v 'fc 6 


crrio2 #i 


IFe/HI 


o-[Fe/H] 


N 


Dm 


0-fDrnl 




Tvoe 


(mag) 


(K) 




(dex) 




(dex) 










HD099998 

x LLJ\jy 7770 


K4III 


+ 1.24 


3863 


60.9 


1.79 


0.18 


-0.16 


0.09 


4 


1.186"'" 


0.011 


HD101501 


G8Vvar 


+3.58 


5401 


60.9 


4.60 


0.18 


-0.13 


0.09 


3 


1.054 


0.009 


HD 102224 


KOIII 


+0.98 


4383 


75.0 


2.02 


0.40 


-0.46 


0.15 


2 


1.148 


0.011 


HD 102328 


K3III 


+2.63 


4390 


60.9 


2.09 


0.18 


+0.35 


0.09 


2 


1.155 


0.011 


HD1 03095 


G8Vp 


+4.37 


5025 


60.9 


4.56 


0.18 


-1.36 


0.09 


4 


1.043"'" 


0.002 


HD103877 

X X XV J. v_/ u / / 


AmV 


+5.88 


7341 


117.6 


4.00 


0.21 


+0.40 


0.09 


1 


1.038 


0.009 


HD 104307 


K2III 


+3.68 


4451 


117.6 


2.00 


0.21 


-0.01 


0.10 


2 


1.159 


0.012 


HD 105262 


B9 


+6.75 


8542 


117.6 


1.50 


0.21 


-1.37 


0.10 


2 


1.056 


0.010 


HD107213 


F8Vs 


+5.13 


6298 


60.9 


4.01 


0.18 


+0.36 


0.09 


4 


1.051 f 


0.009 


HD1 10014 


K2III 


+2.01 


4399 


60.9 


1.47 


0.18 


+0.10 


0.09 


3 


1.151 


0.011 


HD1 1 1631 


M0.5V 


+4.88 


3785 


60.9 


4.75 


0.21 


+0.10 


0.09 


4 


1.109"'" 


0.011 


HD 113092 


K2III 


+2.11 


4283 


60.9 


1.95 


0.18 


-0.37 


0.09 


3 


1.148 


0.011 


HD 113285 


M8III 


-1.06 


2485 


117.6 


0.00 


0.18 






5 


1.250 


0.012 


HD1 14038 

l lis _L i TU J (J 


K1III 


+2.72 


4530 


60.9 


2.71 


0.18 


-0.04 


0.09 


4 


1.132 t 


0.011 


HD 114095 


G5 


+5.88 


4650 


60.9 


2.40 


0.18 


-0.70 


0.09 


2 


1.108 


0.011 


HD1 14330 


AVI 


+4.30 


9509 


117.6 


3.80 


0.21 


-0.01 


0.10 


2 


1.057 


0.011 


HD 114946 


G6V 


+3.11 


5171 


60.9 


3.64 


0.18 


+0.13 


0.09 


3 


1.093 


0.010 


HD 114961* 


M7III 


+ 1.50 


3012 


112.8 


0.00 


0.18 


-0.81 


0.09 


8 


1.295"'" 


0.010 


HD117176 


G5V 


+3.50 


5525 


60.9 


3.81 


0.18 


-0.10 


0.09 


2 


1.067 


0.009 


HD1 17876 


G8III 


+3.87 


4782 


75.0 


2.25 


0.40 


-0.50 


0.15 


3 


1.106 


0.012 


HD 119228 


M2III 


+0.33 


3600 


60.9 


1.60 


0.18 


+0.30 


0.09 


8 


1.203 1 


0.011 


HD 119667 


K5 


+4.18 


3700 


117.6 


1.00 


0.21 


-0.35 


0.10 


4 


1.219"'" 


0.011 


HD120933 


K5III 


-0.00 


3820 


60.9 


1.52 


0.18 


+0.50 


0.09 


8 


1.226 f 


0.012 


HD121130 


M3TTT 


-0.24 


3672 


117.6 


1.25 


0.21 


-0.24 


0.10 


2 


1.216 


0.012 


HD121299 


K2III 


+2.85 


4710 


60.9 


2.64 


0.18 


-0.03 


0.09 


4 


1.125"'" 


0.011 


HD122563 


F8IV 


+3.73 


4566 


75.0 


1.12 


0.40 


-2.63 


0.15 


4 


1.053"'" 


0.011 


HD 122956 


G6IV-Vw 


+5.89 


4635 


75.0 


1.49 


0.40 


-1.75 


0.15 


2 


1.084 


0.010 


HD 123299 


AOIII 


+3.63 


9894 


60.9 


2.90 


0.21 


+0.12 


0.10 


5 


1.048 


0.011 


HD1 23657 


M4III 


-0.23 


3450 


60.9 


0.85 


0.21 


+0.00 


0.09 


4 


1.231"'" 


0.013 


HD124186 


K4III 


+3.59 


4347 


60.9 


2.10 


0.18 


+0.24 


0.09 


3 


1.144 


0.011 


HD 124850 


F7IV 


+2.80 


6116 


60.9 


3.87 


0.18 


-0.11 


0.09 


3 


1.051 


0.009 


HD 124897 


K2IIIp 


-2.91 


4361 


75.0 


1.93 


0.40 


-0.53 


0.15 


2 


1.137 


0.013 


HD 126327 


M7.5III 


+ 1.74 


2819 


60.9 


0.00 


0.18 


-0.58 


0.09 


4 


1.282"'" 


0.013 


HD 126681 

x iiy J. \jvju J. 


G3V 


+7.63 


5536 


60.9 


4.65 


0.18 


-1.25 


0.09 


2 


1.068 


0.010 


HD 126778 


KOIII 


+5.84 


4847 


60.9 


2.34 


0.21 


-0.62 


0.09 


2 


1.120 


0.011 


HD 127243 


G3IV 


+3.15 


4978 


75.0 


3.20 


0.40 


-0.59 


0.15 


3 


1.091 


0.012 


HD 130694 


K4III 


+ 1.10 


4040 


60.9 


1.85 


0.18 


-0.34 


0.09 


3 


1.173 


0.011 


HD 130705 


K4II-III 


+3.95 


4336 


60.9 


2.10 


0.18 


+0.41 


0.09 


4 


1.151 f 


0.011 


HD 13 1430 


K2/K3III 


+2.20 


4190 


60.9 


2.18 


0.18 


+0.04 


0.09 


4 


1.159 f 


0.011 


HD131918 


K4III 


+2.09 


3970 


60.9 


1.49 


0.18 


+0.22 


0.09 


3 


1.173 


0.012 


HD 132345 


K3III-IVp 


+3.26 


4374 


60.9 


1.60 


0.18 


+0.23 


0.09 


3 


1.155 


0.013 


HD 134063 


G5III 


+5.47 


4885 


60.9 


2.34 


0.21 


-0.69 


0.09 


4 


1.109 f 


0.011 


HD 135722 


G8III 


+ 1.22 


4847 


75.0 


2.56 


0.40 


-0.44 


0.15 


2 


1.086 


0.011 


HD 136726 


K4III 


+ 1.93 


4120 


60.9 


2.03 


0.18 


+0.07 


0.09 


4 


1.190 + 


0.014 


HD 137471 


Ml III 


+ 1.03 


3422 


60.9 


1.10 


0.18 


+0.07 


0.10 


3 


1.194 


0.010 


HD 137704 

x xxv I. *j i i \j r 


K4III 


+2.11 


4095 


60.9 


1.97 


0.18 


-0.27 


0.09 


4 


1.1 66 1 


0.011 


HD 137759 


K2III 


+0.77 


4498 


60.9 


2.38 


0.18 


+0.05 


0.09 


3 


1.125 


0.012 


HD 137909 


FOp 


+3.45 


8541 


117.6 


4.25 


0.21 


+0.83 


0.10 


3 


1.047 


0.009 


HD1 38481 

X lly J. «J 0~0 J. 


K5III 


+ 1.21 


3890 


60.9 


1.64 


0.18 


+0.20 


0.09 


4 


1.166"'" 


0.018 


HD1 39641 


G7 5IIIb 


+3.10 


5030 


60.9 


3.22 


0.18 


-0.55 


0.09 


2 


1.081 


0.011 


HD140160 


A0p...V 


+5.20 


9164 


117.6 


3.30 


0.21 


-0.25 


0.10 


2 


1.050 


0.011 


HD141527 


G01ab:pe 


+4.56 


6816 


60.9 


0.48 


0.18 


-0.50 


0.09 


4 


1.042 


0.014 


HD141714 


G3.5III 


+2.66 


5230 


60.9 


3.02 


0.18 


-0.29 


0.09 


5 


1.075 


0.012 


HD 145675 


K0V 


+4.71 


5264 


75.0 


4.66 


0.40 


+0.34 


0.15 


4 


1.080"'" 


0.009 


HD 146051 


M0.5III 


-1.17 


3793 


117.6 


1.40 


0.21 


+0.32 


0.10 


2 


1.189 


0.011 


HD 147923 


M... 


+3.46 


3600 


117.6 


0.80 


0.21 


-0.19 


0.10 


4 


1.221"'" 


0.011 


HD 148783 


M6III 


+0.29 


3279 


112.8 


0.20 


0.21 


-0.06 


0.09 


4 


1.249 f 


0.012 


HD148897 


G8pn 


+ 1.96 


4284 


117.6 


1.15 


0.21 


-0.75 


0.09 


2 


1.116 


0.011 



Table A.l. continued. 



Name 


Spectral 


K 


/ off 

' en 


u L 1 err j 


l02 2 


(xN02 2\ 


IFe/HI 


o-[Fe/H] 


N 


Dm 


0-fDrnl 




Tvne 


(mag) 


(K) 




(dex) 




(dex) 










HD 149009 


K5III 


+2.04 


3910 


60.9 


1.60 


0.18 


+0.30 


0.09 


3 


1.209 


0.012 


X 11_V i i y \J\J _L 


KOV 


+3.91 


5168 


60.9 


4.63 


0.18 


+0.04 


0.09 


4 


1.075"'" 


0.009 


HD150012 


F5III-IV 


+5.26 


6505 


60.9 


3.90 


0.18 


+0.16 


0.09 


3 


1.058 


0.009 


HD150680 


F9IV 


+ 1.28 


5672 


60.9 


3.74 


0.18 


+0.01 


0.09 


2 


1.055 


0.009 


HD151203 


M31IIab 


+0.89 


3640 


60.9 


0.70 


0.18 


-0.10 


0.10 


3 


1.218 


0.011 


HD1 54733 

X XXV _L JT / *J 


K3III 


+2.48 


4279 


60.9 


2.10 


0.18 


+0.00 


0.09 


4 


1.165' 


0.011 


HD1 55763 


B6III 


+3.60 


13397 


117.6 


4.24 


0.18 


-0.95 


0.10 


2 


1.049 


0.010 


HD1 56014 

X I LV 1J\JU1~ 


M51h-ll 

1VU1U XX 


-1.99 


3161 


112.8 


0.00 


0.21 


+0.00 


0.09 


2 


1.267 


0.012 


HD1 56026 


K5V 


+3.86 


4541 


60.9 


4.54 


0.18 


-0.37 


0.09 


2 


1.106 


0.009 


HD1 56283 


K311var 


-0.02 


4460 


60.9 


2.33 


0.18 


+0.18 


0.09 


3 


1.176 


0.011 


HD1 57910 

X IL/ J. «J / 7 1U 


G5III 


+4.39 


5137 


60.9 


1.83 


0.18 


-0.32 


0.09 


3 


1.077 


0.011 


HD1 60933 

X XXV J. \J\J J 


F9V 


+4.82 


5684 


75.0 


3.90 


0.40 


-0.32 


0.15 


2 


1.061 


0.009 


HD161096 

X XXV J. \J J. V v/ 


K2III 


+0.43 


4543 


60.9 


2.16 


0.18 


+0.08 


0.09 


2 


1.138 


0.016 


HD161797 

X XXV 1U1 / v / 


G5IV 


+ 1.51 


5411 


75.0 


3.87 


0.40 


+0.16 


0.15 


3 


1.067 


0.010 


HD161817 


A2VI(HB) 


+6.29 


7759 


60.9 


2.95 


0.18 


-0.95 


0.09 


2 


1.046 


0.010 


HD163990 


M6Svar 


+0.19 


3365 


117.6 


0.70 


0.21 


+0.01 


0.10 


3 


1.247 


0.011 


HD163993 


G81II 


+ 1.73 


5028 


60.9 


2.70 


0.18 


+0.03 


0.09 


3 


1.087 


0.012 


HD 16405 8 


K5III 


-1.16 


3930 


60.9 


1.26 


0.18 


-0.05 


0.09 


4 


1.1 89 1 


0.012 


HD164136 


F2II 


+2.77 


6799 


117.6 


2.63 


0.21 


-0.30 


0.10 


3 


1.044 


0.011 


HD 164349 


K0.5IIb 


+ 1.93 


4446 


60.9 


1.50 


0.18 


+0.39 


0.09 


4 


1.157+ 


0.011 


HD 166208 


G8III... 


+2.93 


4919 


75.0 


2.52 


0.40 


+0.08 


0.15 


3 


1.075 


0.011 


HD 167768 


G3III 


+3.89 


5235 


60.9 


1.61 


0.21 


-0.68 


0.09 


4 


1.085"'" 


0.011 


HD168322 


G8.5IIIb 


+3.93 


4793 


60.9 


2.00 


0.18 


-0.40 


0.09 


3 


1.111 


0.012 


HD1 68720 

X XXV X v/U / \J 


M1III 


+0.74 


3810 


60.9 


1.10 


0.18 


+0.00 


0.10 


4 


1.21 1 1 


0.012 


HD1 68723 

X XXV X v/U / i-^*J 


KOIII-IV 


+ 1.05 


4859 


75.0 


3.13 


0.40 


-0.19 


0.15 


2 


1.099 


0.011 


HD173819 


KOIbpvar 


+2.15 


4421 


117.6 


0.00 


0.21 


-0.88 


0.10 


3 


1.121 


0.012 


HD174638 


B7Ve.... 


+3.19 


12136 


60.9 


2.50 


0.18 


+0.43 


0.09 


3 


1.041 


0.011 


HD175865 

X XXV 1 / *J KJ \J^J 


M5III 


-1.83 


3520 


60.9 


0.50 


0.18 


+0.14 


0.10 


4 


1.250"'" 


0.012 


HD181096 


F6IV: 


+6.47 


6276 


75.0 


4.09 


0.40 


-0.26 


0.15 


2 


1.047 


0.009 


HD182835 


F2Ib 


+4.01 


7350 


721.4 


2.15 


0.32 


+0.09 


0.29 


2 


1.050 


0.010 


HD 184499 


GOV 


+5.07 


5738 


100.0 


4.02 


0.50 


-0.66 


0.30 


4 


1.044 f 


0.009 


HD 184786 


M4.5III 


+0.74 


3467 


117.6 


0.60 


0.21 


-0.04 


0.10 


4 


1.234"'" 


0.012 


HD 185 144 


KOV 


+2.90 


5260 


75.0 


4.55 


0.40 


-0.24 


0.15 


4 


1.061 1 


0.009 


HD187216 


R... 


+6.02 


3500 


117.6 


0.40 


0.21 


-2.48 


0.10 


4 


1.138 t 


0.010 


HD 187921 


KOvar 


+3.80 


6000 


117.6 


1.00 


0.21 


+0.28 


0.10 


3 


1.046 


0.012 


HD188119 


G8III 


+ 1.73 


4915 


75.0 


2.61 


0.40 


-0.32 


0.15 


3 


1.098 


0.012 


HD191277 


K3III 


+2.72 


4459 


60.9 


2.71 


0.18 


+0.30 


0.09 


4 


1.131 + 


0.011 


HD195593 


F5Iab 


+3.65 


6700 


721.4 


1.95 


0.18 


+0.12 


0.09 


3 


1.046 


0.012 


HD 199799 

X XXV 177 / Z> Z> 


Mil 


+ 1.33 


3400 


117.6 


0.30 


0.21 


-0.24 


0.10 


4 


1.241 + 


0.012 


HD202447 


G0III+... 


+2.34 


6087 


60.9 


3.24 


0.18 


+0.09 


0.09 


5 


1.085 


0.010 


HD209369 

X 11VZ.U7 Juy 


F5V 


+3.96 


6217 


60.9 


3.85 


0.18 


-0.26 


0.09 


2 


1.051 


0.009 


HD2 16228 


K0III 


+ 1.27 


4768 


75.0 


2.49 


0.40 


+0.01 


0.15 


2 


1.118 


0.011 


HD217382 


K4III 


+ 1.49 


4035 


60.9 


1.24 


0.18 


-0.25 


0.09 


2 


1.182 


0.011 


HD223047* 


G5Ib 


+2.47 


4990 


117.6 


1.50 


0.21 


+0.18 


0.10 


2 


1.140 


0.010 


HD232078 


K4-5III 


+4.19 


4008 


60.9 


0.30 


0.18 


-1.73 


0.09 


4 


1.139 t 


0.011 



Table A.2. Cluster stars used in the fitting function procedure. AGB stars are maked with a ★ symbol. 



Name 


T 

-* eff 


rr\ Tiff 1 
^ L 1 eit J 


1o£? 2 
Aw e> o 




IFe/HI 

L-i v^/ x xj 


<x[Fe/H] 


Dm 


crFDrnl 




(K) 




(dex) 






(dex) 








Lillerl-6 

1 < 1 1 1 V 1 1 Yt 


3612 


127.0 


+0.05 


0. 


1 1 

1 1 


-0.61 


0.09 


1.250 


0.047 


Lillerl-7 


3671 


96.0 


+0.07 


0. 


13 


-0.61 


0.09 


1.178 


0.043 


Lillerl-157* 


4011 


121.0 


-0.11 


0. 


1 1 

1 1 


-0.61 


0.09 


1.248 


0.047 


Lillerl-158 


3671 


96.0 


-0.02 


0. 


1 1 

1 1 


-0.61 


0.09 


1.245 


0.047 


Lillerl-162 

1 ; 1 1 1 V . ] 1 1 Yt _ 


3627 


119.0 


+0.29 


0. 


1 1 


-0.61 


0.09 


1.214 


0.045 


Lillerl-166 


3973 


118.0 


+0.56 


0. 


13 


-0.61 


0.09 


1.150 


0.041 


Lillerl -299 


3150 


134.0 


-0.07 


0. 


1 1 


-0.61 


0.09 


1.255 


0.047 


M69-1 


3830 


106.0 


+0.04 


0. 


10 


-0.78 


0.03 


1.215 


0.023 


M69-II-37 


3716 


97.0 


+0.09 


0. 


1 ^ 
12 


-0.78 


0.03 


1.178 


0.022 


M69-I-40 


3917 


113.0 


+0.25 


0. 


1 ^ 
12 


-0.78 


0.03 


1.166 


0.022 


M69-2 


3864 


109.0 


+0.30 


0. 


1 1 


-0.78 


0.03 


1.215 


0.023 


M69-3 


3864 


109.0 


+0.30 


0. 


1 1 


-0.78 


0.03 


1.212 


0.023 


M69-4 


3899 


112.0 


+0.41 


0. 


1 ^ 
12 


-0.78 


0.03 


1.219 


0.024 


M71-29 


3641 


108.0 


+0.09 


0. 


12 


-0.84 


0.06 


1.203 


0.017 


M71-30 


3992 


120.0 


+0.68 


0. 


13 


-0.84 


0.06 


1.185 


0.016 


M71-B 


3764 


100.0 


+0.25 


0. 


12 


-0.84 


0.06 


1.188 


0.016 


M71-46 


4011 


121.0 


+0.80 


0. 


1 1 


-0.84 


0.06 


1.195 


0.016 


M71-A4 


4153 


134.0 


+0.79 


0. 


13 


-0.84 


0.06 


1.178 


0.016 


M71-1=H 


3796 


103.0 


+0.86 


0. 


1 O 


-0.84 


0.06 


1.175 


0.016 


M7 1-2=1 


4565 


173.0 


+ 1.08 


0. 


1 1 


-0.84 


0.06 


1.149 


0.015 


M71-3=113 


3954 


116.0 


+0.83 


0. 


13 


-0.84 


0.06 


1.192 


0.016 


M7 1-4=45 


4011 


121.0 


+0.87 


0. 


13 


-0.84 


0.06 


1.171 


0.016 


M71-5=64 


4153 


134.0 


+ 1.42 


0. 


10 


-0.84 


0.06 


1.159 


0.016 


M7 1-6=66 

1VX / X \J \J\J 


4649 


182.0 


+ 1.39 


0. 


10 


-0.84 


0.06 


1.137 


0.015 


M71-8=21 


4458 


162.0 


+ 1.45 


0. 


13 


-0.84 


0.06 


1.138 


0.015 


NGC0104-A02 

1 ' V J V^ V./ X V 1^ 1 \ V./ Z— 


3533 


128.0 


+0.08 


0. 


13 


-0.78 


0.02 


1.214 


0.035 


NGC0104-W12 


3780 


102.0 


+0.05 


0. 


1 l\ 

10 


-0.78 


0.02 


1.225 


0.036 


NGC0104-A19 

1 1 V.J V^ V/ X U i z\. X J/ 


3554 


128.0 


+0.10 


0. 


12 


-0.78 


0.02 


1.181 


0.034 


NGC0104-V07 

1 1 V.J V^ V./ X V l V V/ / 


3764 


100.0 


+0.21 


0. 


13 


-0.78 


0.02 


1.188 


0.034 


NGC0104-V06 


3864 


109.0 


+0.43 


0. 


12 


-0.78 


0.02 


1.172 


0.033 


NGC0104-L168 

1 < v J V^ V./ X \J^T \—i \-\J\j 


3882 


110.0 


+0.43 


0. 


12 


-0.78 


0.02 


1.215 


0.035 


NGC01 04-5529 


3973 


118.0 


+0.64 


0. 


13 


-0.78 


0.02 


1.205 


0.035 


NGCO 104-2426 


4070 


126.0 


+0.93 


0. 


1 O 

12 


-0.78 


0.02 


1.170 


0.033 


NGCO 104-1 505 


4070 


126.0 


+0.96 


0. 


1 ^ 
12 


-0.78 


0.02 


1.184 


0.034 


NGC01 04-441 8 

1 ' V J v.- V/ X V./ 1^ I I X YJ 


4091 


128.0 


+0.95 


0. 


l l 


-0.78 


0.02 


1.188 


0.034 


NGC0104-1510 


4153 


134.0 


+ 1.01 


0. 


1 A 

14 


-0.78 


0.02 


1.168 


0.033 


NGC01 04-241 6 

1 i v.J V^ V./ X V/ z— X V./ 


4219 


140.0 


+ 1.33 


0. 


15 


-0.78 


0.02 


1.177 


0.034 


NGCO 104-6408 


4383 


155.0 


+ 1.44 


0. 


14 


-0.78 


0.02 


1.163 


0.033 


NGC0288-A96 


4070 


126.0 


+0.50 


0. 


10 


-1.14 


0.03 


1.144 


0.016 


NGC0288-A78 


4132 


132.0 


+0.67 


0. 


1 1 


-1.14 


0.03 


1.148 


0.016 


NGC0288-C20 

1 i v.J V ^ V./ Z— v.J U V ^ ^— V./ 


4153 


134.0 


+0.71 


0. 


1 1 


-1.14 


0.03 


1.118 


0.015 


NGC0288-A77 


4219 


140.0 


+0.83 


0. 


1 A 

10 


-1.14 


0.03 


1.135 


0.016 


NGC0288-A245 

i i v j v ^ v./ z— v.J u y \ z— 7 


4433 


160.0 


+ 1.10 


0. 


12 


-1.14 


0.03 


1.124 


0.015 


NGC0362-III1 1 

i 1 VJ VyV*.' Ui- XXX X X 


4011 


121.0 


+0.37 


0. 


1 A 

10 


-1.09 


0.03 


1.162 


0.035 


NGC0362-IV100 


4091 


128.0 


+0.50 


0. 


1 A 

10 


-1.09 


0.03 


1.145 


0.034 


NGC0362-TTT63 

1 1 VJ V--- V./ *J YJ 1 1 1 V J . ' 


4031 


123.0 


+0.52 


0. 


1 A 

10 


-1.09 


0.03 


1.169 


0.036 


NGC0362-TTT44 

i i v.J v^ v./ ^ V/ ill i r 


4091 


128.0 


+0.66 


0. 


1 A 

10 


-1.09 


0.03 


1.119 


0.033 


NGC0362-III70 

1 1 V_J V^ V/ YJ XXX / vy 


4310 


148.0 


+0.75 


0. 


1 A 

10 


-1.09 


0.03 


1.088 


0.031 


NGC5927-100* 


3847 


107.0 


+0.01 


0. 


1 1 


-0.64 


0.01 


1.246 


0.046 


NGC5927-799 


3847 


107.0 


+0.38 


0. 


1 O 

12 


-0.64 


0.01 


1.243 


0.045 


NGC5927-627 


3864 


109.0 


+0.51 


0. 


l l 


-0.64 


0.01 


1.203 


0.043 


NGC5927-532 


3992 


120.0 


+0.80 


0. 


12 


-0.64 


0.01 


1.167 


0.041 


NGC5927-622 


4175 


136.0 


+ 1.00 


0. 


15 


-0.64 


0.01 


1.152 


0.040 


NGC5927-536 


4310 


148.0 


+ 1.27 


0. 


15 


-0.64 


0.01 


1.191 


0.042 


NGC6388-1 


3954 


116.0 


-0.34 


0. 


13 


-0.74 


0.18 


1.238 


0.038 


NGC6388-3 


3899 


112.0 


-0.05 


n 


1 1 
1 1 


-0.74 


0.18 


1.230 


0.038 


NGC6388-4 


3701 


96.0 


+0.09 


0. 


13 


-0.74 


0.18 


1.204 


0.037 


NGC6440-KF-1 


3847 


107.0 


+0.41 


0. 


13 


-0.62 


0.10 


1.179 


0.018 


NGC6440-KF-2 


3813 


105.0 


+0.46 


0. 


11 


-0.62 


0.10 


1.210 


0.018 



Table A.2. continued. 



Name 








o-[\ogg] 


[Fe/H] 


<r[Fe/H] 


Deo 


o-[D co ] 




(K) 




(dex) 






(dex) 










NGC6440-KF-3 


3716 


97.0 


+0.47 


o 


12 


-0.62 


0.10 


I 


1 96 


0.018 


NGC6440-KF-4 


3686 


94.0 


+0.50 




12 


-0.62 


0.10 


I 


193 


0.018 


NGC6440-KF-5 


3864 


109.0 


+0.53 


n 


1 1 


-0.62 


0.10 


I 


239 


0.019 


NGC6440-KF-6 


3747 


99.0 


+0.57 




12 


-0.62 


0.10 


I 


1 88 


0.018 


NGC6440-KF-8 


3747 


99.0 


+0.67 


o 


XX 


-0.62 


0.10 


I 


1 96 


0.018 


NGC6440-KF-7 


3796 


103.0 


+0.67 


n 


12 


-0.62 


0.10 


I 


.201 


0.018 


NGC6528-7 


3864 


109.0 


+0.30 


n 


1 1 


-0.61 


0.08 


I 


.220 


0.029 


NGC6528-11 


3732 


98.0 


+0.05 





1 1 


-0.61 


0.08 


I 


.242 


0.030 


NGC6528-22* 


3813 


105.0 


+0.37 


n 


12 


-0.61 


0.08 


I 


275 


0.031 


NGC6528-6* 


3917 


113.0 


+0.41 




12 


-0.61 


0.08 


I 


.248 


0.030 


NGC6553-20 


3780 


102.0 


-0.16 




1 ? 


-0.60 


0.04 


1 




0.034 


NGC6553-19 


3551 


138.0 


-0.02 


o 


10 


-0.60 


0.04 


I 


.227 


0.035 


NGC6553-25 


3747 


99.0 


+0.14 


o 


14 


-0.60 


0.04 


1 


.229 


0.036 


NGC6553-16* 


3917 


113.0 


+0.27 


n 


XX 


-0.60 


0.04 


1 


288 


0.038 


NGC6553-26* 


3847 


107.0 


+0.32 


o 


1 3 


-0.60 


0.04 


I 


.281 


0.038 


NGC6553-14* 


3813 


105.0 


+0.38 




12 


-0.60 


0.04 


1 


265 


0.037 


NGC6553-2 


3656 


105.0 


+0.42 


n 


12 


-0.60 


0.04 


1 


192 


0.034 


NGC6624-KF-1 


3917 


113.0 


+0.40 


o 


1 1 


-0.70 


0.03 


I 


230 


0.034 


NGC6624-KF-2 


3917 


113.0 


+0.53 


o 


1 3 


-0.70 


0.03 


I 


1 60 


0.030 


NGC6624-KF-3 


3917 


113.0 


+0.54 


n 


1 3 


-0.70 


0.03 


I 


.162 


0.030 


NGC6624-KF-4 


3899 


112.0 


+0.82 




14 


-0.70 


0.03 


1 


.218 


0.033 


NGC6624-KF-5 


4070 


126.0 


+0.96 






-0.70 


0.03 


I 


1 97 

. ±y i 


0.032 


NGC6712-LM5* 


4111 


130.0 


+0.55 


n 


10 


-0.94 


0.03 


I 


.223 


0.024 


NGC6712-LC01 


4132 


132.0 


+0.58 




1 1 


-0.94 


0.03 


I 


.211 


0.023 


NGC6712-LC03 


4219 


140.0 


+0.61 




10 


-0.94 


0.03 


1 


1 88 

. 1 OO 


0.023 


NGC6712-LM8 


4111 


130.0 


+0.71 


o 


1 1 


-0.94 


0.03 


I 


1 57 


0.022 


NGC6712-LM10 


4264 


144.0 


+0.72 




1 


-0.94 


0.03 


I 


194 


0.023 


NGC6712-B66 


4196 


138.0 


+0.85 




1 ? 

1 z. 


-0.94 


0.03 


1 


. i yu 


0.023 


Terzan2-1 


4241 


142.0 


+0.40 


o 


10 


-0.65 


0.14 


I 


1 

.1 JO 


0.034 


Terzan2-2 


3973 


118.0 


+0.41 




1 1 


-0.65 


0.14 


I 


205 


0.037 


Terzan2-3* 


3899 


112.0 


+0.43 


o 


12 


-0.65 


0.14 


I 


?60 


0.040 


Terzan2-4 


4175 


136.0 


+0.44 


0. 


11 


-0.65 


0.14 


1 


.215 


0.037 


Terzan2-5 


3936 


115.0 


+0.56 


0. 


12 


-0.65 


0.14 


1 


.200 


0.037 


Terzan2-7 


4111 


130.0 


+0.61 


0. 


10 


-0.65 


0.14 


1 


.199 


0.037 


Terzan2-8 


4132 


132.0 


+0.62 


0. 


10 


-0.65 


0.14 


1 


.151 


0.034 
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